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Abstract 


This work is the continuation of à longstanding project which aims at developping 
a fixed-wing unmanned air vehicule. This UAV has à high-mounted wing, a pushing 
propeller and its designated payload is a stereoscopic camera fixed on the front of the 
fuselage, surrounded by a fixed plexiglas globe. 


Based on the four previous years of developpement, two colleagues and I have modified 
and developped the design further. This year, a real-scale wing was manufacturerd in order 
to conduct mechanical tests on it. Furthermore, a new, aerodynamic-efficient half-scale 
fuselage was developped to study the aerodynamics of the new design in the wind tunnel 
of the University of Liège. 


This report is articulated about four main chapters. The first develops the changes 
applied to the modified fuselage, the constraints we faced in conceiving it and experi- 
mental data gathered from the wind tunnel. Second, an aerodynamic study is conducted 
on the whole half-scale model, that is including a wing and an empennage. Alas, due to a 
technical problem with a 6-axes sensor of the wind tunnel, no equilibrium study has been 
conducted ; only 3-axes mesures were available. Third, an in-depth study is conducted 
on the propulsive system. The motor was also tested in the wind tunnel in order to 
obtain the quadratic thrust curve, and hence the performances at take off and in flight. A 
design change to a twin-engine configuration is also introduced. Last, the flight control 
surfaces are sized following given flight constraints, and the servomotors acting on those 
surfaces are selected based on the couple computed at the hinge of the mobile surface. 


As such, this work tries to address some hurdles raised in the previous reports, parti- 
cularly on the aerodynamic and performance side. 
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Chapitre 1 


Introduction 


Ce travail de fin d’études fait partie d’un projet qui consiste à concevoir et développer 
un drone. Les caractéristiques et critères de design ont constamment évolué au cours des 
quatre dernières années. Les bases de la construction de ce projet seront résumées dans 
les lignes qui suivent, afin que le lecteur soit averti du jargon et des concepts utilisés par 
après. 


1.1 Historique 


Le projet de développement, conception et construction d’un drone de prises de vues 
en 3D a été initié par le Département d’Aérospatiale et Mécanique de l’Université de Liège, 
en collaboration avec le service du Prof. Jacques G. Verly du département Montéfiore de 
la même université. 


Année académique 2009 - 2010 Durant cette première année, un groupe de quatre 
étudiants à réalisé un prédimensionnement — dont les traits principaux sont représentés 
à la figure 1.1(a) — qui permet de répondre à la mission proposée, à savoir prendre des 
photographies en 3D de bâtiments ou autres lieux d'intérêt. Ces photographies sont prises 
en inscrivant le drone dans une trajectoire circulaire centrée sur la cible. 


Par ailleurs, une étude détaillée de la charge utile est menée : design d’un mécanisme 
de support de la caméra 3D et test de cette dernière, étude stéréoscopique et définition des 
paramètres permettant d'obtenir la meilleure image. Subséquemment, a pu être réalisé un 
choix des moteurs et réducteurs orientant la caméra selon 2 axes. Ainsi, une géométrie 
particulière sera proposée, ce qui permet d'établir l'enveloppe de vol en fonction de la 
vitesse maximale et de l’altitude maximale. 

À également été réalisée durant cette première année, une étude aérodynamique de 
base qui a permis de dimensionner les surfaces portantes et l’empennage. Les profils choisis 
pour les ailes seront respectivement pour l’aile et pour les surfaces de l’empennage, un 
NACA 2412 et des NACA 0012. Ce choix fut motivé par une volonté de se baser sur des 
profils bien connus, largement utilisés et documentés. 


L'étude inclut également le design d’un système propulsif. Il à été recherché un moteur 
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simple à mettre en place, puissant et minimisant les problèmes dynamiques tels que les 
vibrations, néfastes pour la qualité des prises de vue. Le choïx fut ainsi porté vers un 
moteur électrique issu du monde de l’aéromodélisme : un moteur brushless alimenté par 
des batteries de type LiPo. 


Une partie expérimentale fut finalement menée afin de comparer, sur un modèle à 
l'échelle 1:2, les propriétés aérodynamiques mesurées et les résultats numériques. 


Année académique 2010 - 2011 Une grande partie du travail réalisé en 2010 a été 
de résoudre les problèmes aérodynamiques mis en évidence lors des essais réalisés l’an- 
née précédente. Les surfaces de commandes sont ainsi dimensionnées pour répondre aux 
contraintes de design et améliorer la stabilité latérale. De plus, la section des poutres liant 
l’'empennage à l’aile (et donc au fuselage) a été changée (d’une section carrée à une section 
circulaire) afin d'obtenir de meilleures propriétés mécaniques. Les changements apportés 
sont illustrés à la figure 1.1(b). 


Une deuxième partie a permis de déterminer les exigences structurales à supporter par 
le drone. Un certain nombre de modifications ont été apportées à la maquette initiale de 
2010, par exemple afin de permettre un réglage de l’angle de calage de l’empennage. 


L'ensemble de ces modifications ayant entraîné une augmentation notoire de la masse 
totale, il fallut compenser ce poids par une plus grande portance. Dès lors, il fut décidé 
de rallonger l’aile (passant ainsi de 3 à 3,3 mètres) ainsi que d'augmenter la vitesse de 
vol (de 55 à près de 70 km/h). Des essais ont été réalisés sur le moteur afin de vérifier 
que la poussée fournie soit suffisante, ce qui n’a pas été le cas : une analyse a montré que 
la traînée seule ne pouvait être vaincue par la poussée, sans doute suite à une mauvaise 
circulation de l’air autour du fuselage et réduisant ainsi l’efficacité de l’hélice. 


Année académique 2011 - 2012 Sur base des observations formulées en 2011, il est 
décidé de modifier la forme du fuselage, dans l’espoir de diminuer la traînée — et d’aug- 
menter le flux d’air visible par le moteur. Toutefois, les modifications apportées n’ont pas 
sensiblement réduit la traînée mesurée en soufHerie. Une visualisation de l’écoulement en 
soufflerie, par l’application de tufts, a cependant montré que l’écoulement sur la maquette 
à l’échelle 1:2 semblait être meilleur. 


Par ailleurs, la conception de l’aile —- désormais d’un seul tenant fixée au fuselage, 
plutôt que deux demi-ailes — ainsi que sa structure interne et les éléments raccordant les 
poutres de l’'empennage au fuselage, a été développée en détail (choix des tôles et des 
colles, montage, etc). Les matériaux utilisés sont également modifiés : le composite en 
fibres de carbone remplace le bois et la structure interne est intégralement en aluminium. 


Des tests en soufflerie sont également menés sur la maquette du nouveau fuselage 
pour permettre une étude des stabilités longitudinale et latérale. Une étude détaillée du 
placement du centre de gravité et son influence sur les grandeurs aérodynamiques est 
entreprise sur base du nouveau design. Ce dernier est illustré à la figure 1.1(c). 
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1.2 Situation en 2012-2013 


Dans l’ensemble des travaux précédents, un ensemble de perspectives est proposé et 
un certain nombre de questions restent à solutionner. Les tâches à réaliser cette année ont 
été réparties entre trois étudiants. 


Norman Ernst s’est penché sur la construction de l’aile et, auparavant, sur l’usinage 
des moules positifs et négatifs ainsi que sur les matériaux, colles etc. permettant d’obtenir 
une aile en fibre de carbone à l’échelle 1:1. Des essais mécaniques ont été effectués sur 
des échantillons de composites, et un test de traction effectué sur l’aile a montré que 
la construction de l’aile, effectuée par Norman et l’auteur avec l’aide des techniciens du 
laboratoire MDF, était prometteuse. L’aile construite n’est cependant plus utilisable, car 
le test a endommagé la structure. De plus amples détails sont disponibles en consultant 
la référence [2]. 


Vladimiro Capizzi s’est quant à lui intéressé à la structure (interne) du fuselage ainsi 
qu’au dimensionnement du train d’atterrissage. Un exosquelette est ainsi proposé, com- 
prenant un support pour le mécanisme entraînant la caméra, sur base des contraintes mé- 
caniques calculées par une méthode éléments finis dans SAMCEF . Le train d’atterrissage 
est dimensionné à partir d’un bilan énergétique. De plus amples détails sont disponibles 
en consultant la référence [6]. 


Ce travail, basé sur les caractéristiques reprises à la table 1.1, reprend les avancées 
effectuées au niveau du fuselage, et étudie l’effet du changement de géométrie sur les 
propriétés aérodynamiques tant du fuselage (chap. 2) que du drone complet (chap. 3). Il 
étudie également les performances du drone (chap. 4). Des essais en soufflerie sur le moteur 
ont ainsi permis de caractériser la poussée du moteur et d’en dériver les caractéristiques au 
décollage et en vol de l’appareil. Enfin, on trouvera dans ce travail un prédimensionnement 
des surfaces de commande ainsi que des servomoteurs (chap. 5). 


Aiïle | Envergure 3,3 m 
Surface 1,419 m° 
Corde 43 cm 
Allongement 7,7 
Profil NACA 2412 
Surface de l’aileron 0,08 m? 
Empennage horizontal | Envergure 1m 
Surface 0,2 m°? 
Corde 20 cm 
Profil NACA 0012 
Surface de l’élévateur | 510 cm? 
Empennage vertical | Envergure 37 cm 
Surface 740 cm? 
Corde 20 cm 
Profil NACA 0012 
Surface du gouvernail | 204 cm? 
Dimensions globales | Longueur du fuselage | 1,15 m 
Longueur totale 2,13 m 
Hauteur & 97 cm 
Masse estimée 14 kg 


TABLE 1.1 - Caractéristiques du drone. Cas de référence dans ce travail. 
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1.3 Enjeux de ce travail 


D'abord, un des défis était la réalisation d’un modèle à l’échelle 1:1 d’une maquette 
de fuselage permettant de tester les performances du moteur en conditions quasi réelles. 
Il s’est avéré que l’étape de fabrication a été très lente, notamment suite aux contraintes 
matérielles (machines hors services) et humaines (techniciens très demandés). Ainsi, seul 
un modèle à l’échelle 1:2 à été réalisé. Les contraintes de fabrication nous ont imposé une 
modification du fuselage, ce qui a nécessité de nouveaux tests en soufflerie pour valider le 
modèle. 


Par ailleurs, une étude théorique, numérique et expérimentale a été effectuée sur le 
drone complet. À cette fin, la théorie des plaques minces a été utilisée en parallèle avec 
le code TORNADO et des essais en soufflerie afin de comparer les différentes grandeurs 
entre elles. Seront expliquées en détail la modélisation par éléments finis pour le calcul 
numérique, ainsi que les corrections effectuées pour les mesures en soufHerie. 


Ensuite, un des objectifs du présent travail est l’étude des performances du moteur. 
Le moteur testé a été acquis les années précédentes et ne répondait pas aux attentes, 
notamment suite à la traînée trop importante du fuselage et une puissance insuffisante. 
Des essais ont donc été réalisés afin de pousser le moteur au maximum de ses capacités, et 
établir les courbes de poussée quadratiques relatives à ces essais. À partir de ces mesures, 
il à été possible de s'intéresser aux performances requises tant au décollage que pour 
certains cas de vol. Une introduction au cas bimoteur et ses effets inhérents au décollage 
est également proposée. 


Enfin, un autre objectif de ce travail était le dimensionnement des surfaces mobiles de 
commande, ainsi que la mise en œuvre pratique de la transmission du couple aux surfaces 
de contrôle. Une étude tant de type analytique que numérique est réalisée afin de comparer 
le dimensionnement des surfaces, alors qu’une méthode analytique/empirique est utilisée 
pour le calcul des couples aux charnières. Sur base des résultats précédents, une ébauche 
de système de contrôle est proposée et analysée d’un point de vue pratique. 
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3300 


(c) Drone 2012 (d) Drone 2013 


FIGURE 1.1 - Évolution du drone au cours des années. Infographie reprise de [2]. 


Chapitre 2 


Fuselage 


Ce deuxième chapitre est consacré à l’étude de la géométrie et du comportement aéro- 
dynamique du fuselage. Une première partie rappelle brièvement les étapes qui ont mené 
au design proposé dans les travaux de 2012. Ensuite, il sera explicité un ensemble des 
contraintes qui s’appliquent sur le fuselage ; ces contraintes vont mener à apporter des 
modifications et à développer une ultime esquisse de fuselage. Une étude expérimentale 
de la maquette du drone sera menée en soufflerie afin de comparer le comportement du 
nouveau fuselage à celui de 2012. 


+ La construction du fuselage à l’échelle 1:2 est reprise dans le travail [2], le design 
et le développement de la structure interne du fuselage à l’échelle 1:1 sont étudiés dans le 
travail [6]. 


2.1 Evolution de la géométrie, d’année en année 


Cet historique succinct se base sur les travaux réalisés les années précédentes. L’en- 
semble de ces travaux sont repris dans la bibliographie. 


La première année, un prédimensionnement avait mené à la réalisation d’une première 
maquette à l’échelle 1:2 dont les formes et la construction étaient rudimentaires. La figure 
2.1 montre la maquette durant un essai en soufferie. 


Une maquette à l’échelle 1:1 fut ensuite réalisée : le fuselage est constitué de panneaux 
plans en bois et la sphère comportant la caméra est remplacée par un hémisphère en frigo- 
lite suspendu sous le fuselage. À l’arrière du fuselage, on trouve le moteur et ses batteries. 
L'intérieur du fuselage comprend l’électronique de commande ainsi que deux cloisons de 
séparation, faisant office de frames. Après des études approfondies en soufflerie, il s’est 
avéré qu'un tel fuselage, couplé avec le moteur choisi!, ne permettait pas d'obtenir une 
poussée suffisante pour voler et, même, de décoller. Certaines modifications furent appor- 
tées, notamment sur le dimensionnement des gouvernes et surfaces de commandes pour 
améliorer le comportement aérodynamique, mais pas de remise en question fondamen- 


1. On trouvera les propriétés de ce moteur à l’annexe D.2.1 


CHAPITRE 2. FUSELAGE A 


FIGURE 2.1 — Illustration de la maquette 2010 (échelle 1:2). 


tale du fuselage. La figure 2.2 montre la maquette à l’échelle 1:1 avec son équipement 
embarqué, toujours dans la soufHerie. 


FIGURE 2.2 -— Illustration de la maquette 2011 (échelle 1:1). 


En vue de solutionner la poussée insuffisante, ou plutôt la traînée trop importante, 
une modification fondamentale du fuselage fut envisagée. Les parois furent arrondies et la 
sphère, comportant la caméra 3D, remontée pour se dissimuler dans le fuselage. Cepen- 
dant, contre toute attente, les essais en soufflerie n’ont pas montré une diminution nette 
de la traînée induite par le fuselage. Il est dès lors supposé que le design même du fuselage 
pose problème et doit être revu. L'intégration de la sphère, par exemple, devrait être la 
plus lisse possible et perturber l’écoulement au minimum, ce qui n’est pas le cas avec 
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la géométrie de base acceptée jusqu'ici. La dernière version du fuselage, une maquette à 
l'échelle 1:2, est illustrée à la figure 2.3. 


FIGURE 2.3 — Illustration de la maquette 2012 (échelle 1:2). 


2.2 Contraintes sur le fuselage 


Sur base de l’historique développé ci-avant, il apparaît qu’une optimisation du fuse- 
lage s'impose. Les modifications dépendent de quatre paramètres principaux. Le premier 
paramètre à rencontrer est la capacité machine à pouvoir usiner le fuselage. Le second 
concerne le mécanisme qui permet de manœuvrer la caméra en vol. Le troisième et ultime 
paramètre est la conservation de la qualité de prise de vue, soit essentiellement de garder 
une fenêtre de vue la plus large possible. 


2.2.1 Considérations techniques 


Le drone est réalisé entièrement dans le Laboratoire de Méthode de Fabrication de 
l’Université de Liège, dirigé par le Prof. Debongnies. La conception des pièces est donc 
largement tributaire de l'équipement disponible et fonctionnel du laboratoire. Hélas, du- 
rant cette année académique, nous n’avons pas pu bénéficier du concours de la fraiseuse 
multi-axes Huron, ce qui aura pour effet de fortement ralentir les phases d'usinage des 
moules et maquettes, et dès lors de limiter notre temps et les phases d'assemblage du 
drone. 


Le fuselage proposé en 2012 est certes très esthétique, mais il ne répond hélas que 
partiellement aux critères d'amélioration émis en 2011. De plus, suite à l’indisponibilité 
du centre d'usinage Huron, il est apparu que les moules (positifs) de ce fuselage ne savaient 
pas être usinés en l’état — sinon de manière très coûteuse en temps, en ressources et en 
argent. 
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Dès lors, il a été décidé d’apporter des modifications au fuselage, améliorations qui 
répondent à la fois au besoin d'améliorer l’aérodynamique selon des procédés de bonne 
pratique, et à la capacité machine que nous utiliserons. Dans un souci de facilité, nous 
utiliserons principalement une fraiseuse 3-axes à commande numérique de la société Mi- 
kron ?. 


Le principal inconvénient du fuselage proposé en 2012 est, en ce qui concerne l’uti- 
lisation de la Mikron, une course de la table insuffisante pour un usinage efficace et 
économique du fuselage, plus particulièrement de la partie avant, c’est-à-dire la sphère 
contenant la caméra stéréoscopique. Une solution très simple et pragmatique a donc été 
de remonter la sphère par rapport à l’axe du fuselage. 


2.2.2 Considérations aérodynamiques 


Le fuselage doit être optimisé pour réduire au plus la traînée et/ou améliorer l’écou- 
lement de l’air autour du moteur et augmenter ainsi le rendement de propulsion. Cette 
amélioration se fait selon des principes généraux et du bon sens. Les résultats issus des 
essais en soufHerie sont détaillés à la section 2.4. 


2.2.3 Considérations mécaniques 


La caméra stéréoscopique est orientée dans une direction à l’aide d’un mécanisme 
composé de moteurs et réducteurs. Dès lors, le fuselage proposé doit permettre d’enfermer 
l’ensemble de ce mécanisme. À cet effet, nous avons profité de la venue d’un représentant de 
la société Minimotor, de Herentals, qui nous à proposé le meilleur compromis en moteurs 
et réducteurs répondant aux exigences exposées dans [7]. 


A l’aide des datasheets de ces éléments, nous avons pu déterminer que la dimension 
initiale de la sphère, soit un diamètre de 30 cm, convenait parfaitement. Elle ne sera donc 
pas modifiée. Par ailleurs, l’existence de moteurs compacts à permis de fortement limiter 
la nécessité d’avoir un espace important au-dessus de la sphère. 


L'ensemble des éléments précédents nous confirment l’intérêt de remonter la sphère, 
mais il reste encore un point, et non des moindres, à éclaircir. 


2.2.4 Considérations optiques 


La première géométrie de fuselage permettait d'obtenir une large fenêtre d'observation, 
de l’ordre de 360 degrés mais avec les trains d’atterrissage comme obstacles. Cette fenêtre 
a été réduite durant l’année 2012 afin de permettre l'intégration de la sphère au fuselage. 
Pour le nouveau fuselage que nous proposons, il apparaît que la fenêtre se réduit à un angle 
de 94 degrés, comme l’illustre la figure 2.4. Il s’agit certes d’une dégradation notable de 
la qualité du design initial, mais en l’état actuel et compte tenu des contraintes énoncées 
précédemment, il nous semble que la nouvelle géométrie est un bon compromis entre 
aérodynamique et fabrication. 


2. Les caractéristiques de la fraiseuse CNC Mikron sont reprises à l’annexe D.1.1 
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FIGURE 2.4 — Disposition de la caméra dans la demi-sphère ; fenêtre d’observation pour 
la Tachyon 3D. 


2.3 Élaboration d’un nouveau fuselage 


Sur base des contraintes exposées ci-dessus, on a dessiné un nouveau fuselage sur le 
logiciel CATIA. La géométrie du nouveau fuselage est illustrée à la figure 2.5, tandis 
qu’une photographie représentative est proposée à la figure 2.6. Les valeurs numériques 
sont des mesures de longueur en millimètres. 


La modification principale apportée dans ce travail par rapport au fuselage de 2012 
concerne, on le voit, la position de la sphère frontale. En effet, toute la partie arrière, soit 
les 700 mm en partant de l’extrémité de la queue, ainsi que la position de l’aile ont été 
conservés par rapport au design précédent. Cependant, puisque le fuselage est légèrement 
plus court, le centre de gravité va se déplacer vers l’arrière. Pour le conserver à la moitié de 
la corde de l’aile et à mi-hauteur du fuselage, il faudra jouer sur la position des batteries 
et des systèmes de commande. 


2.4 Comparaison des fuselages 
Essais en soufflerie 


Nous allons à présent mener une étude comparative des derniers fuselages développés. 
Le tableau 2.1 reprend les surfaces frontales des trois fuselages (échelle 1:1), illustrées à 
la figure 2.7. Les comparaisons seront effectuées sur bases des données de [8]. 


Les tests seront effectués pour deux vitesses de vent : 10 et 20 m/s (consigne). Les 
mesures seront effectuées avec un mât instrumenté * pour calculer les trois composantes de 


3. La calibration du mât est effectuée avant chaque nouvelle expérience en soufflerie. La calibration 
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300 mm 


64 mm |250 mm 


115 m 


(a) Vue du haut. 


(b) Vue latérale. 


FIGURE 2.5 — Illustration de la nouvelle géométrie du fuselage adoptée pour la suite du 
projet. Cotes en millimètres. 


FIGURE 2.6 - Maquette du fuselage à l’échelle 1:2 


Fuselage | Surface frontale (m?) 


2011 0,12 
2012 0,096 
2013 0,071 


TABLE 2.1 — Surface frontale des fuselages proposés les dernières années. 


sera détaillée à la section 3.3.2. Par ailleurs, il faudra veiller à supprimer la traînée induite par le mât. 
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Fuselage 2011 


Fuselage 2012 


Fuselage 2013 


FIGURE 2.7 — Surface frontale des fuselages proposés les dernières années. 


forces appliquées sur le fuselage. Pour chaque configuration d’angle d’attaque et de vent, 
trois mesures sont effectuées. Lors du traitement numérique, une moyenne sera effectuée 
entre ces trois mesures et seules les grandeurs moyennes de C', et C seront considérées. 


Les grandeurs de référence pour le calcul de ces deux coefficients, à partir des données 
expérimentales mesurées sur le fuselage à l’échelle 1:2, sont : 

— Température : 18,5° ; 

— Pression atmosphérique : 982 hPa ; 

— Masse volumique de l'air : 1,211 kg/m° : 

— Surface : 0,355 m°? 


2.4.1 Traînée 


L'évolution du coefficient de traînée Ch en fonction de l’angle d'attaque du fuselage 
a est illustrée à la figure 2.8. Pour rappel, le coefficient Ch est défini comme par : 


D 


a 
+ Ce 


(2.1) 
où D est la force de portance mesurée et corrigée, p la masse volumique de l’air, U la 
vitesse du vent mesurée au tube Pitot et S,, la surface de l’aile. 


On constate que la courbe du nouveau fuselage suit assez fidèlement celle du fuselage 
de 2011 et reste généralement inférieure à celle du fuselage de 2012 (sauf autour de à — 
4°). Notons cependant que le support permettant le maintien de la maquette au mât de 
la soufferie, illustré à la figure 2.6, crée une traînée parasite au niveau de l’aile. Il faut 
également se rendre compte que la variation du coefficient C7 est assez faible et devient 
ainsi sensible aux erreurs de mesure. 


En extrapolant les résultats de Ch à 19,9 m/s pour le fuselage à taille réelle, on peut 
montrer que la traînée est comprise entre 7 et 11 N. Par rapport au fuselage des années 
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FIGURE 2.8 — Coefficient de traînée du fuselage en fonction de l’angle d’attaque. 


précédentes, il apparaît donc que l'optimisation de la traînée est en bonne voie (sauf 
autour de a — 4°). 


2.4.2 Portance 


Par analogie avec le paragraphe précédent, l’évolution du coefficient de portance C% 
est illustrée à la figure 2.9. Contrairement aux fuselages précédents, on constate que la 
nouvelle géométrie crée une certaine quantité de la portance, et ce quel que soit l’angle 
d'attaque étudié. Cette génération de portance est certes de faible amplitude, mais sur le 
fuselage à taille réelle et pour une vitesse de 19,9 m/s, elle peut varier de 2 à 11 N selon 
l’angle d'attaque considéré. 


0.05 


FIGURE 2.9 - Coefficient de portance du fuselage en fonction de l’angle d’attaque. 
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La figure montre également que la portance augmente globalement avec l’angle d’at- 
taque, ce qui est avantageux pour la portance totale générée par le drone au décollage — 
mais n’en tirons pas de conclusions trop hâtives. 


2.4.3 Écoulement de l’air autour du corps 


Afin de pouvoir visualiser sommairement l’écoulement autour du nouveau fuselage, 
nous avons placé à plusieurs endroits stratégiques des morceaux de coton dont une ex- 
trémité peut se mouvoir librement dans l’air. Ces tufts seront placés sous la sphère ainsi 
que sur le cône de queue. La figure 2.10 montre l’écoulement autour des fuselages 2012 et 
2013. Dans les deux cas, l’écoulement semble bien attaché au fuselage, sauf juste derrière 
la sphère où l’écoulement a un comportement tourbillonnaire. 


La zone sensible, devant le moteur, semble épargnée par les tourbillons. Il s’agit d’un 
point crucial pour obtenir le meilleur rendement de propulsion de l’hélice. 


(a) Fuselage 2012 


AT 


(b) Fuselage 2013 


FIGURE 2.10 — Illustration de l’écoulement par des tufts. 
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2.5 Synthèse de ce chapitre 


À partir du design proposé en 2012, nous avons modifié la position de la sphère par 
rapport à l’axe du fuselage. Une maquette à l’échelle 1:2 a été réalisée et testée en soufflerie. 
Les résultats montrent que les objectifs à améliorer sont grosso modo atteints : la traînée 
est en moyenne diminuée par rapport aux fuselages des années 2011 et 2012, tandis que 
l'écoulement sur le cône de queue reste convenable. En bonus, le nouveau fuselage crée 
une quantité limitée de portance à tout angle d’attaque testé. 


Chapitre 3 


Aérodynamique 


Ce chapitre est consacré à la modélisation de l’aérodynamisme du drone. Dans un 
premier temps, on recherche des ordres de grandeur pour les coefficients de portance 
Cx et de moment C; à partir de la théorie des profils minces. Ensuite, une simulation 
plus détaillée prenant en compte la géométrie et les propriétés des ailes est effectuée à 
l’aide du code TORNADO. Enfin, nous pourrons comparer ces résultats avec des valeurs 
expérimentales obtenues en soufHerie. 


3.1 Théorie des profils minces 


La théorie des profils minces suppose que l’écoulement soit bidimensionnel (le profil 
d’allongement infini), stationnaire et incompressible. 


Le rapport d’allongement de l’aile (Aspect Ratio) dans le plan x-y est À = 7,674 et 
son rapport de conicité est de 1. En connaissant le profil NACA de l’aile principale, on 
peut déterminer l’angle d’attaque à portance nulle de l’aile, à savoir ao. On peut trouver 
dans [9] (pp. 321-322) que ap = —2, 1° pour un profil NACA 2412. 


La pente du coefficient de portance pour des profils minces est ao = dat À 27, comme 
expliqué à l’annexe A. 


3.1.1 Coefficient de portance 


Le coefficient de portance total alors peut s’obtenir par les équations suivantes : 


2T 
E =1+ ————— .l 
Fr À(1+7) (1) 
__ 0,995a 
daCL = Era (3.2) 
CL — dCL(a = Q) (3.3) 


où À est le rapport d’allongement et 7 le rapport de conicité. 


L’angle de calage de l’aile sur le fuselage, noté 4,,, est de trois degrés. Dès lors, l’angle 
d'attaque total de l’aile a, = à +1, est de trois degrés lorsque le fuselage est à incidence 
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nulle. 


3.1.2 Coefficient de traînée 


Le coefficient de traînée peut se calculer à partir d’une approximation de la polaire de 
la traînée. L’équation de la polaire s'exprime par : 


CRC (3.4) 


où Cn, est la traînée indépendante de la portance et Ÿ est le facteur d’Oswald. Ces 
deux valeurs peuvent être approximées à partir de la table 3.1. Dans le cas du drone, on 
supposera que l'efficacité du fuselage se situe entre celle d’un avion monomoteur à train 
fixe et d’un avion agricole, soit Cn, — 0,045 et Ÿ — 0,75. 


Aircraft Type (Cr e 
High-subsonic jet 0.014-0.020 0.75-0.85 
Large turboprop 0.018-0.024 0.80-0.85 
Twin-engine piston 0.022-0.028 0.75-0.80 
aircraft 

Single-engine piston 0.020-0.030 0.75-0.80 
aircraft with fixed gear 

Single-engine piston 0.025-0.040 0.65-0.75 
aircraft with retractable 

gear 

Agricultural aircraft 0.060 0.65-0.75 
without spray system 

Agricultural aircraft with  0.070-0.080 0.65-0.75 


spray system 


TABLE 3.1 - Valeurs communément admises pour Ch, et Ÿ. Source : [5] 


3.1.3 Coefficient de tangage 


On peut également calculer une valeur de la pente de la droite représentant C,, en 
fonction de l’angle d’attaque. Partant de la figure 3.1, on peut écrire l’expression générale 
du coefficient de tangage pour un empennage dénué de surfaces mobiles : 


Cm = bo + (ao + a a)(h — ho) (3.5) 


où &o, 41 et bd sont des coefficients connus, À est la distance entre le centre de gravité et 
le bord d’attaque de l’aile eths est la distance entre le centre aérodynamique et le centre 
de gravité. On obtient alors immédiatement que : 


daCm = a(h — ho) < 0 (3.6) 


où a; est la pente de portance de l’empennage en fonction de l’angle d’attaque «, soit, à 
nouveau en se basant sur un profil mince, & 27. La grandeur d,C,, doit être négative pour 
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he 
FIGURE 3.1 — Forces et moments appliqués sur l’appareil en tangage. Figure issue de [5]. 


avoir une configuration stable. En effet, après une perturbation de l’angle d'incidence a, 
l’accroissement de portance doit être contrebalancé par le poids appliqué au centre de 
gravité. Cette condition impose ainsi que hQ soit plus petit que h, d’où le signe négatif. 


Multipliant ce coefficient par le bras de levier À — h9 et par l’angle d’attaque, on peut 
facilement déterminer le coefficient de tangage, supposant que les coefficients ao et bo 
soient nuls. 


3.2 Modélisation sur TORNADO 


Pour poursuivre notre étude de l’aérodynamique du drone, nous faisons appel à un 
code MATLAB opensource, TORNADO. Le lecteur intéressé par ce code MATLAB peut 
consulter la référence [10]. 


TORNADO permet de calculer les grandeurs aérodynamiques (ainsi que leurs dérivées) 
des ailes et des empennages, et prend en compte l’angle de déflexion des surfaces de 
commande. Dès lors, les effets liés au fuselage ne peuvent pas être modélisés. Il en résultera 
une sous-estimation non négligeable de la traînée totale du modèle, mais la portance 
calculée devrait être assez proche de la portance réelle, puisque le fuselage ne participe 
presque pas à la création de portance. 


3.2.1 Dimensionnement du drone 


Peu de modifications ont été apportées par rapport à la géométrie proposée par [11]. 
Seule, on l’a vu, la partie frontale s'étendant du sommet de la sphère au bord d’attaque 
de l’aile est modifiée. Aïnsi, le modèle aérodynamique des ailes est identique à celui utilisé 
par [8]. La géométrie des ailes qui va être modélisée dans TORNADO est rappelée à la 
figure 3.2. 
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200 mm 


500 mm 


430 mm 


FIGURE 3.2 — Positionnement et longueur des ailes et surfaces de commandes. 


3.2.2 Modélisation des surfaces portantes 


Quatre modèles/configurations vont être étudiés sur TORNADO. 
1. l’aile principale; 

2. l’élévateur ; 

3. le gouvernail ; 

4. l’assemblage des surfaces formant le drone. 


La modélisation TORNADO utilise des éléments finis pour le calcul des différentes 
grandeurs. Il nous faut dès lors rechercher la meilleure répartition de panneaux: celles qui 
minimise le temps de calcul et l’erreur de calcul commise. Cette distribution sera choisie 
linéaire: chaque partition (sous-ensemble de l’aile modélisée, permettant de définir les 
zones avec et les zones sans surface mobile) est donc représentée par un ensemble donné de 
panneaux identiques. Nous pourrons alors obtenir un ensemble de grandeurs relativement 
fiables nous permettant de connaître plus précisément le comportement aérodynamique 
de chaque surface prise individuellement. 
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Aile principale 


L’aile du drone est symétrique et de distribution de corde constante. Elle a une enver- 
gure de 3,3 mètres et une corde de 43 centimètres. Par la théorie des profils minces, on 
sait que le centre aérodynamique de l’aile est situé au voisinage du quart de corde, soit à 
10,75 cm du bord d’attaque. Le point de référence utilisé pour le calcul des moments est 
superposé au centre de gravité de l’aile. 


Pour cette partie de l’étude aérodynamique, la dimension des ailerons sera fixée à 
55 cm d'envergure et 40% de la corde. L’aile va alors être divisée en deux partitions 
symétriques !. Le première partition a une envergure de 1,1 m, tandis que la seconde fait 
55 cm et définit l’envergure de l’aileron. Ces deux partitions définissent une aile sans 
dièdre ni flèche, dont le rapport de conicité est unitaire. Le calage ne sera pas modifié et 
sera de 3°. 


Le tableau 3.2 et la figure 3.3 reprennent les configurations de panneaux qui vont être 
testées. La figure 3.4 illustre l’effet du nombre de panneaux sur les grandeurs aérodyna- 
miques lors d’une étude paramétrique de la déflexion (non symétrique) des gouvernes, 
pour un angle Ô variant de -45° à 45° par pas de 5°. Cette figure montre que, pour mi- 
nimiser l'erreur sur C7, Cp et Cm, il est préférable d'utiliser un maillage de type 5 par 
rapport à l'erreur commise et aux ressources informatiques requises. 


Panneaux sur partition intérieure Panneaux sur partition extérieure 
Config. | Corde | Enverg. | Corde gouverne | Corde | Enverg. | Corde gouverne 
Type 1 2 4 0 2 2 1 
Type 2 3 5 (0 3 2 2 
Type 3 4 8 0 4 3 2 
Type 4 5 9 (0) 5 3 2 
Type 5 6 10 0 6 4 2 
Type 6 8 12 0 8 6) 3 


TABLE 3.2 - Types de configurations utilisées pour modéliser l’aile principale. 


FIGURE 3.3 — Types de configurations utilisées pour modéliser l'aile. 


1. Les grandeurs qui suivent doivent donc être multipliées par 2 pour représenter fidèlement le modèle 
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FIGURE 3.4 - Grandeurs aérodynamiques en fonction du nombre de panneaux et de l’angle 
de déflexion des ailerons, 6. 


Empennage vertical 


L’empennage proposé pour le drone comprend deux dérives avec leur gouverne de 
direction, reliées au reste du fuselage par deux poutres de section circulaire s’enchâssant 
dans l’aile portante. Le profil sélectionné pour réaliser la dérive est un NACA0012, soit un 
profil symétrique de rapport longueur/épaisseur de 12%. Les envergures précédemment 
établies sont de 37 cm pour l’aile et la gouverne. De plus, les cordes respectives sont de 
20 cm et 6 cm. Ces valeurs seront à revoir en fonction des résultats complets fournis par 


l'analyse TORNADO. 


Le tableau 3.3 et la figure 3.5 reprennent les configurations de panneaux qui vont être 
testées. La figure 3.6 illustre l’effet du nombre de panneaux sur les grandeurs aérody- 
namiques lors d’une étude paramétrique de la déflexion des gouvernes, pour un angle Ô 
variant de -45° à 45° par pas de 5°. Cette figure montre que, pour minimiser l’erreur sur 
Ch, il est préférable d'utiliser un maillage de type 2 par rapport à l'erreur commise et 
aux ressources informatiques requises. 
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Panneaux sur seule partition 
Config. | Corde | Enverg. | Corde gouverne 


Type 1 1 2 1 
Type 2 2 4 2 
Type 3 4 6 3 


TABLE 3.3 - Types de configurations utilisées pour modéliser l’empennage vertical. 


Pa | 
Type 1 Type 2 Type 3 


FIGURE 3.5 —- Types de configurations utilisées pour modéliser l’empennage vertical. 
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FIGURE 3.6 —- Grandeurs aérodynamiques en fonction du nombre de panneaux et de l’angle 
de déflexion de la gouverne, 0. 
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Empennage horizontal 


L’empennage horizontal est situé entre les deux dérives. Le profil ainsi que les cordes 
de l’aile et de la gouverne de profondeur sont identiques à celles de l’empennage vertical. 
Cependant, on va créer deux partitions ? : la première comprend l’élévateur et a une 
envergure de 42,5 cm tandis que la seconde a une envergure de 7,5 cm. 


Le tableau 3.4 et la figure 3.7 reprennent les configurations de panneaux qui vont être 
testées. La figure 3.8 illustre l’effet du nombre de panneaux sur les grandeurs aérody- 
namiques lors d’une étude paramétrique de la déflexion de l’élévateur, pour un angle Ô 
variant de -45° à 45° par pas de 5°. Cette figure montre que, pour minimiser l’erreur sur 
Cp et Cm, il est préférable d'utiliser un maillage de type 2 par rapport à l’erreur commise 
et aux ressources informatiques requises. 


Panneaux sur partition extérieure 
Corde | Enverg. | Corde gouverne 


Panneaux sur partition intérieure 
Config. | Corde | Enverg. | Corde gouverne 


Type 1 2 2 1 2 1 (0) 
Type 2 3 3 2 3 2 (0) 
Type 3 4 4 3 4 3 (0 


TABLE 3.4 — Types de configurations utilisées pour modéliser l’empennage horizontal. 
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FIGURE 3.7 — Types de configurations utilisées pour modéliser l’empennage horizontal. 


2. cf. supra 
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FIGURE 3.8 — Grandeurs aérodynamiques en fonction du nombre de panneaux et de l’angle 
de déflexion de l’élévateur, 6. 


3.2.3 Modélisation du drone 


Le comportement des différentes surfaces est certes intéressant, mais notre ambition est 
de faire voler non une aile mais un drone. À cet effet, l’ensemble des surfaces portantes est 
modélisé selon la géométrie finale. TORNADO ne permet pas, dans son code de calcul, 
la prise en compte du fuselage et de ses propriétés. Aïnsi, les grandeurs relatives à la 
traînée sont largement sous-estimées, au contraire des grandeurs relatives à la portance 
étant donné que le fuselage contribue très peu à la portance totale. 


Le tableau 3.5 résume les différents paramètres utilisés pour modéliser le drone dans 
TORNADO, tandis que la figure 3.9 illustre le modèle à l’échelle avec son maillage. Les 
résultats issus de cette modélisation seront exploités dans la suite du travail. 


Entité | Envergure | Corde | Enverg. gouv. | Corde gouv. | Maillage 
Aile 330 cm 43 cm | 50 cm 15,5 cm Type 5 
Empennage horiz. | 100 cm 20 cm | 85 cm 6 cm Type 2 
Empennage vert. | 37 cm 20 cm | 37 cm 6 cm Type 2 


TABLE 3.5 - Modèle TORNADO utilisé pour l’étude aérodynamique du drone. 


Dans tous les calculs effectués, le centre de gravité ainsi que le point de référence sont 
placés, comme conseillé par [8] : 

— horizontalement, à 50 % de la corde de l’aïle ; 

— verticalement, à 50 % de la hauteur du fuselage (sans la sphère). 
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— LL 


Vue babord Vue frontale 


Vue en plan Vue 1sométrique 


FIGURE 3.9 —- Modèle éléments finis du drone. 


3.3 Essais en soufflerie 


Le capteur multi-axes NANO25, utilisé pour mesurer les différents moments aérody- 
namiques, s’est avéré être non-fonctionnel au moment de mener les essais en soufferie. 
Puisque l’ensemble de la maquette à été dimensionné pour fonctionner avec ce capteur 
(masse maximale, couples appliqués, etc.) et au vu des contraintes de temps et d’argent 
disponible, il ne sera malheureusement pas possible de faire des essais de stabilité dans 
ce travail. Dès lors, les concepts de stabilité latérale et longitudinale seront basés sur les 
essais menés en 2012. 


3.3.1 Équilibrage du modèle réduit 


Avant toute expérience en soufferie, il s’avère nécessaire d’équilibrer au maximum la 
maquette complète, à savoir le fuselage, l’aile et l’empennage, afin que l’ensemble des 
couples au centre de gravité soient nuls, ceci afin d'éviter la saturation des capteurs au 
repos. À cet effet, il a été prévu, lors de la conception du support permettant d’attacher 
la maquette au mât de la soufflerie, un système de masse mobile permettant de jouer 
finement sur la position du centre de masse que l’on veut situer au niveau du mât. 


Un petit montage pratique, illustré à la figure 3.10, nous a permis de calculer le couple 
à appliquer pour contrebalancer le poids de l’empennage. Pour cela, nous avons attaché 
une poutre au support de l’empennage, et placé à l’extrémité de la poutre un seau que 
l’on remplit d’eau. Connaissant le point d'application ainsi que le poids du seau et de son 
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contenu, on peut calculer la position et la masse à appliquer en tout point situé sur l’axe 
longitudinal. 


CG 


FIGURE 3.10 —- Montage permettant de déterminer le lest à appliquer au nez pour équili- 
brer la maquette 1:2. 


La position longitudinale à l’avant de la maquette étant réduite, on essaye d’appliquer 
l’ensemble des efforts à une distance maximale de 25 cm en amont du centre de gravité. 
Par un bilan d'équilibre, on détermine qu’une masse de 2,7 kg appliquée à une distance 
d'environ 24 cm du mât permettait de maintenir la maquette en équilibre sur son pied. 


On conçoit un support permettant de relier le fuselage, l’aile et le fuselage au mât de 
la soufferie. Le support comprend un tube vertical, au-dessus duquel on soude un plateau 
dans lequel vient se loger le barre maintenant l’aile et le fuselage, et en dessous duquel 
on fixe une rondelle permettant de fixer le tube au capteur NANOZ25. Un tube horizontal 
sera soudé sur le tube vertical. Dans ce tube horizontal, on glissera un autre tube auquel 
le lest est attaché. La position précise du lest peut ainsi être réglée par une vis qui vient 
serrer le tube interne. Afin de pouvoir positionner le fuselage de manière rigide, on colle 
deux rondelles en bois, sur lesquelles la peau sera vissée. Le support, puis le support avec 
la masse mobile sont illustrés à la figure 3.11. 


FIGURE 3.11 — Support permettant la fixation de la maquette complète à l’échelle 1:2 au 
mât de la soufflerie. 


3.3.2 Corrections appliquées à l’écoulement 


Contrairement à l’étude du fuselage, qui est de dimension réduite et a une surface 
frontale plus petite, l’étude de l’écoulement autour du drone complet ne peut se faire sans 
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l’ajout de corrections. En effet, l’envergure de l’aile de la maquette à l’échelle 1:2 du drone, 
soit 1,65 m., est telle que les parois de la souflerie, dont la largeur est de 2 m, ne sont plus 
à une distance très grande du centre de la section. Dès lors, les corrections porteront sur 
l'élimination des effets induits par la soufflerie sur le modèle. 


Calibration du mât La première opération à effectuer avant toute mesure est la 
calibration du mât instrumenté de la soufflerie. Cette opération s'effectue en testant la 
force mesurée par le mât en disposant sur celui-ci des poids de masse connue. La calibration 
du mât dans les directions x et z est reprise en annexe, à la figure B.1. 


Élimination de la traînée du mât Les forces aérodynamiques mesurées, particuliè- 
rement la traînée, sont surestimées par rapport à la réalité puisque l'effet du mât sur 
l'écoulement n’est pas négligeable. Il faut dès lors soigneusement mesurer la traînée du 
mât lors d’un essai à vide pour différentes vitesses de vent et en extraire le coefficient Ch 
correspondant. 


Effet de flottabilité La couche limite qui se développe sur la paroi de la soufflerie tend 
à rétrécir la surface de travail, ce qui a pour effet d’augmenter la traînée mesurée et donc 
une chute de la pression statique. On peut calculer le gradient de pression longitudinal 
par : 


ce (3.7) 


où À est la hauteur de la section de travail et k est un facteur compris entre 0,016 et 
0,04 qui doit être mesuré pour une soufflerie spécifique. Ici, on considère que sa valeur 
est minimale, soit 0,016. La correction de l’effet de flottabilité dépend de ce gradient de 
pression ainsi que du volume du corps VZ de manière telle que : 


VD De LDVB (3.8) 


On calcule le volume du fuselage à l’aide de CATTA, tandis que les ailes sont approximées 
par des parallélépipèdes. Le volume du corps total pour le fuselage à taille réelle est alors 
de 0,276 + 0,074 — 0,35 m°, soit pour la maquette à l'échelle 1:2, 0,044 m°. 


3.4 Résultats et comparaisons 


Afin de valider nos résultats théoriques et numériques, nous les confrontons à des essais 
en soufHerie. 


3.4.1 Stabilité longitudinale 


Le drone sera en position (longitudinalement et statiquement) stable s’il revient spon- 
tanément à sa position initiale après une perturbation de l’angle d’attaque qui l’écarte de 
sa position d'équilibre, soit la configuration pour laquelle il vole de manière stable sans 
intervention sur les commandes. En l’absence des mesures de moment, seules les trois 
forces aérodynamiques seront traitées dans cette section. 
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Moment de tangage 


Ce paragraphe à propos du moment de tangage ne semble pas pertinent lorsque l’on 
sait que nous n'avons pas obtenu de résultats nouveaux par rapport à [8] suite à l’indis- 
ponibilité du capteur multiaxes. Toutefois, il nous permet de définir un angle de calage 
de l’empennage horizontal qui trouvera son utilité tant lors des mesures en soufflerie qui 
suivent que lors du dimensionnement des surfaces de commandes et des servomoteurs, 
réalisé au chapitre 5. 


Une comparaison du moment de tangage en fonction de l’angle d'attaque pour diffé- 
rents angles de calage de l’empennage est alors effectuée entre les résultats expérimentaux 
de 2012 et une analyse TORNADO. Ces résultats sont représentés à la figure 3.12, où les 
traits continus sont les résultats obtenus par TORNADO), les traits discontinus sont les 
résultats (moyennés) de 2012, et la courbe en pointillés est une droite dont la pente est 
donnée par l’équation : 


DaCm = dath — RO) = —2r(h — h0) 


On trouvera le graphique complet de l’essai en soufflerie de 2012 à l’annexe B.2. 


FIGURE 3.12 - Moment de tangage : valeurs théorique, numériques et expérimentales, 
pour plusieurs angles de calage de l’empennage. 


La conclusion issue de [8] propose un angle de calage pour l’empennage de deux degrés. 
Comme le montre la figure 3.12, il ne s’agit pas de la meilleure solution si l’on en croit la 
modélisation numérique. Compte tenu de la précision très relative du calage pouvant être 
réglé en soufflerie, l’angle de calage retenu sera finalement |, = 2°|. Cet angle sera utilisé 
pour les mesures effectuées dans la suite du travail. 


En ce qui concerne la stabilité longitudinale, cette dernière sera rencontrée dès lors 
que dc, Cm < 0. La courbe de C,, en fonction de C7 est représentée à la figure 3.13, 
d’aspect très similaire à la figure 3.12 et pour cause puisque C7, ne dépend que de l’angle 
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d’attaque &, et montre que la stabilité longitudinale est rencontrée tant que le coefficient 
de portance est supérieur à 0, ou encore à un angle d’attaque plus grand que -5°. 


0.2 


—+— Cm 3° num 
— +-Cm3°exp|| 


FIGURE 3.13 — Courbe de C, en fonction de C7. Stabilité longitudinale. 


Portance 


Le coefficient de portance C7 en fonction de l’angle d’attaque du fuselage a est re- 
présenté à la figure 3.14. Cette figure montre que l’évolution est bien linéaire, quelle que 
soit la méthode de calcul. De plus, chaque courbe est relativement proche l’une de l’autre. 
On peut cependant différencier les courbes bleue (plaque plane) et verte (TORNADO) 
des autres (soufferie), les premières tendant à sous-évaluer la portance additionnelle du 
fuselage et les effets 3D sur l’aile que l’on n’a pas pris en compte. 


Un autre facteur important à noter et qui ne se repère que par des essais en soufferie est 
le décrochage. On remarque que pour un angle d’attaque donné (l'angle de décrochage), 
la portance commence à chuter. Dans le cas du drone 2013, cet angle est situé à 6° pour 
un vent de 10 m/s, et à 9° pour un vent de 20 m/s. 


En extrapolant ces résultats au drone à l’échelle 1:1 et en situation réelle, on peut 
obtenir la portance du drone dans plusieurs cas de vol. On peut alors tracer des courbes 
de la portance en fonction de la vitesse, pour plusieurs angles d’attaque, utilisant la 
relation : : 

U 
ie CL Su (3.9) 


avec p — 1,225°. 


De telles courbes sont illustrées à la figure 3.15 et permettent par exemple de vérifier 
que le drone peut voler en croisière. La masse du drone actuel est évaluée à 14 kg, soit 
un poids de 137 N. On détermine alors la vitesse de décrochage : il s’agit de l'intersection 


de la courbe a — 9.2° avec le poids du drone, ce qui correspond à | V, & 12,5 m/s| Ce 
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[= + -0C, num 
a C. exp, U = 10 m/s 
= C, exp, U = 20 m/s 
—+-: C, exp, 2012 


FIGURE 3.14 — Coefficient de portance en fonction de l’angle d’attaque. 


résultat est également vérifié par la relation : 
2W 
Ne 
DCR D 


Traînée 


Le coefficient de traînée C', en fonction de l’angle d’attaque du fuselage a est représenté 
à la figure 3.16. La figure montre un résultat qui peut sembler surprenant, à savoir que le 
CD mesuré à 20 m/s est plus faible que le Cn mesuré à 10 m/s. Il s’agit en fait de l’effet 
de la traînée induite, qui tend à diminuer inversement au carré de la vitesse du vent. 


De plus, on peut voir qu’à l’angle de décrochage, soit a, > 6° ou 9° selon la vitesse du 
vent, la traînée augmente brutalement. 


Enfin, on constate que la traînée calculée de manière numérique, c’est-à-dire issue d’un 
calcul TORNADO, sous-estime largement la traînée mesurée puisque ne prenant pas en 
compte l’effet du fuselage, tandis que la courbe obtenue de manière théorique a été calibrée 
intuitivement avec des paramètres Ch, et Ÿ cohérents, en fonction des caractéristiques du 
drone. 


En extrapolant ces résultats au drone à l’échelle 1:1 et en situation réelle, on peut 
obtenir la traînée du drone dans plusieurs cas de vol. La courbe de la portance en fonction 
de la vitesse, illustrée à la figure 3.17 pour plusieurs angles d'attaque, se trace en utilisant 
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FIGURE 3.15 — Portance totale générée par le drone en fonction de 
(m/s) et de l’angle d'attaque (°). 
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FIGURE 3.16 — Coefficient de traînée en fonction de l’angle d’attaque. 


la relation : 


P 


2 
D = Plus, 


2 


(3.10) 


avec p — 1,225°. À la vitesse de décrochage V, — 12,5 m/s, on aura donc une traînée 
d'environ 13,6 newtons. Cette figure montre bien, si l’on en doutait encore, que la force 
de traînée devient très importante passé l’angle de décrochage. Les résultats obtenus ici 


CHAPITRE 3. AÉRODYNAMIQUE 32 


seront de la plus haute utilité dans l’étude des performances réalisée au chapitre 4. 
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FIGURE 3.17 — Traînée totale générée par le drone en fonction de la vitesse du vent (m/s) 
et de l’angle d’attaque (°). 


Dérapage 


Le dérapage du drone devrait être nul puisque le drone est orienté dans la direction du 
vent. Toutefois, la difficulté à obtenir un réglage fin de l’angle de lacet en soufflerie, ainsi 
que la difficulté à rendre l’élévateur parfaitement plan autour de l’axe du fuselage, va créer 
une force de dérapage, certes de petite amplitude, mais pas nécessairement négligeable, 
comme le montre la figure 3.18 qui représente le coefficient de dérapage Cs en fonction 
de l’angle d’attaque du fuselage a. 


Diagrammes polaires de traînée 


À partir des données précédentes, on réalise deux diagrammes polaires de la portance 
vis-à-vis de la traînée. Le premier diagramme (figure 3.19) est réalisé pour un angle de 
calage de l’'empennage de 3,6° (angle effectif, vu par rapport à l'écoulement), tandis que 
le second (figure 3.20) est réalisé pour un angle de calage de 2,1°. On notera que l’angle 
de downwash n’est pas pris en compte dans les figures proposées. Les polaires sont tracées 
pour deux nombres de Reynolds * : Re — 145270 équivaut à une vitesse de 10 m/s, Re — 
290540 équivaut à une vitesse de 20 m/s. 


Comparant les courbes du graphe 3.19) à celles du graphe 3.20), on constate que les 
résultats obtenus sont similaires quel que soit l’angle de calage 42. Le rebond constaté à la 


Us 
2 


QI 


3. Le nombre de Reynolds de l’écoulement autour d’un profil d’aile s'écrit : Re, = 
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FIGURE 3.18 — Coefficient de dérapage en fonction de l’angle d’attaque. 
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FIGURE 3.19 — Diagramme polaire C7 vs Cn pour deux vitesses avec un calage tn — 3,6°. 


figure (a) résulte de deux mesures effectuées à des angles & presque égaux mais mesurés 
au début et à la fin des essais en soufflerie. 


Dans les deux cas, le décrochage se manifestera pour un C7 proche de l’unité et un 
Ch dix fois moindre. On remarque aussi que la traînée diminue légèrement avec un calage 
plus élevé. Aïnsi, la traînée avec un calage de l’empennage de 2,1° est supérieure à celle 
mesurée en 2012 pour un même calage. Nous supposons qu'il s’agit d’une conséquence 
directe de l'installation en soufflerie, où il fut difficile de régler le calage suite à l’absence 
de gabarit permettant une mesure précise de cet angle. Malgré celà, les valeurs mesurées 
pour un calage de 2,1° sont utilisées pour le calcul des performances : nous nous plaçons 


CHAPITRE 3. AÉRODYNAMIQUE 34 


Polaire, angle de calage (empennage) = 2,1° 


Re = 145270 
Re = 290540 || 
- —- — Essai 2012 


0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 


FIGURE 3.20 — Diagramme polaire C7 vs Ch pour deux vitesses avec un calage 4 — 2,1°. 


ainsi dans le cas défavorable. 


3.4.2 Stabilité latérale 


Par analogie avec l'équilibre longitudinal, le drone sera en position (latéralement) 
stable s’il revient spontanément à sa position initiale après une perturbation des angles 
de lacet et/ou roulis qui l’écarte de sa position d'équilibre. L'étude de la stabilité latérale 
en temps que telle n’est pas très pertinente ici dans la mesure où il n’y a aucun résultat 
expérimental neuf depuis les résultats rapportés dans [8]. Cependant, rappelons que le 
drone sera statiquement stable en lacet et en roulis respectivement si : 


Nous allons tester le comportement de la maquette du nouveau drone en l’orientant 
avec un certain angle de lacet par rapport au vent. On pourra dès lors déjà obtenir un 
aperçu de l’évolution de la portance en fonction de cet angle de lacet. Les autres forces 
aérodynamiques ne sont pas reprises, car le calibrage du mât n’a pas été effectué selon 
l’axe y. 


Portance 


Le coefficient de portance C7, en fonction de l’angle de lacet 5 et en fonction de l’angle 
d'attaque à, est représenté à la figure 3.21. 


Cette figure montre que la portance est relativement stable avec l’angle de lacet pour un 
angle d’attaque donné. Cette stabilité est d'autant plus évidente à grande vitesse de vent. 
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FIGURE 3.21 - Coefficient de portance en fonction des angles d’attaque a et de lacet B. 


Il n’y à d’ailleurs, comme on peut s’y attendre, aucun effet positivant sur le décrochage à 
angle d'incidence trop grand puisque la courbe rouge (pour a — 11°) est systématiquement 
inférieure à la courbe magenta (pour & — 6 ou 9°, angle de décrochage). 


3.5 synthèse de ce chapitre 


Ce chapitre s’est attelé à comparer le comportement aérodynamique du drone complet 
sur base de trois méthodes. La première, la théorie des profils minces, à permis d’obtenir 
des premiers ordres de grandeur. La seconde, une modélisation numérique sur TORNADO), 
a nécessité de définir le modèle éléments finis utilisé pour les calculs afin de minimiser aussi 
bien les erreurs commises que le temps de calcul. La troisième, des essais en soufflerie, 
nous à contraint, d’une part, à effectuer un équilibrage du drone afin de ne pas saturer 
inutilement l’instrumentation ; d'autre part, à considérer des facteurs de correction sur les 
mesures permettant de mieux refléter l’écoulement réel. 


Nous avons ensuite comparé les résultats obtenus par ces trois méthodes. Des résultats 
incomplets obtenus, il ressort que la stabilité longitudinale du drone devrait être assurée 
pour un angle de calage de 2°, tandis que les forces de portance et de traînée, mesurées 
en fonction de l’angle d’attaque, restent proches des essais antérieurs. Par ailleurs, on a 
déterminé que la vitesse de décrochage était de 12,5 m/s. Tout ceci nous laisse penser que 
le comportement du drone à l’échelle 1:1 devrait être acceptable. Les marges de stabilités 
ne pourront être déterminées qu’à partir de mesures de moments en soufflerie. 


Chapitre 4 


Performances 


L'étude aérodynamique présentée au chapitre précédent a permis de déterminer le 
comportement général du drone en vol, mais ne fournit aucune information utile quant 
à la propulsion de l’appareil. Ce chapitre va dès lors s'intéresser à la performance du 
moteur ainsi qu'aux performances du drone lors des phases de décollage et d'atterrissage. 
Nous expliquerons d’abord le montage expérimental du moteur, testerons en soufflerie les 
performances de ce dernier sans fuselage et pour plusieurs cas de charge, étudierons à partir 
de la courbe de poussée ainsi obtenue les performances au décollage, extrapolerons les 
résultats dans le cas d’un drone bimoteur, et enfin discuterons de l’ensemble des résultats. 


4.1 Étude en soufflerie de la poussée du moteur 


4.1.1 Montage expérimental 


Il convient de présenter en quelques lignes le montage effectué en soufferie pour tester 
le moteur. Le montage est illustré à la figure 4.1 et est expliqué ci-après. 


Les expériences qui suivent vont se baser sur l’utilisation d’un moteur électrique bru- 
shless ! de type Graupner Compact 6407. Ce moteur sera alimenté soit par des batteries 
LiPo, soit des batteries 12V. Ces batteries seront telles que la tension d’entrée nominale 
aux bornes du contrôleur du moteur soit de 36-37 V. Le moteur est relié au mât ins- 
trumenté par l’intermédiaire d’un support en bois. Sur ce support sont placés tous les 
éléments qui permettent de faire tourner le moteur, soit le contrôleur, les batteries de 
puissance (LiPo ou 12V), le relais électrique ?, le récepteur (transmission de la consigne à 
partir d’une télécommande) et les batteries alimentant le récepteur (4 piles de 1,2 V). 


Les câbles reliant tous ces composants entre eux seront soigneusement collés ou fixés 
à l’aide de colsons, tout comme la longueur des câbles reliant le moteur aux batteries de 
puissance seront choisis les plus courts possible afin d’éviter un effet Joule trop important. 


On note par ailleurs que l’armature de la soufflerie introduit des perturbations dans 
la transmission radio entre la télécommande et le récepteur. Il sera donc important que 


1. Les caractéristiques techniques de ce moteur sont reprises à l’annexe D.2.1. 
2. Ce composant n’est pas indispensable pour faire tourner le moteur, mais il permet de relier 
distance" le moteur aux batteries de puissance, afin de réaliser des tests en toute sécurité. 


"à 
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FIGURE 4.1 — Schéma du montage expérimental permettant de tester les performances 
du moteur. 


l’antenne soit bien tendue et que la commande ait son antenne à proximité — une solution 
empirique a montré qu'il était préférable de faire passer l’antenne directement dans la 
soufHerie, par un orifice existant. 


Le montage du moteur sur son support est illustré à la figure 4.2, dans le cas d’une 
configuration tirante et alimentation par batteries LiPo. 


FIGURE 4.2 — Illustration du moteur et de son hélice, en configuration tirante, dans la 
soufferie de l’ULg. 
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4.1.2 Comparaison des montages expérimentaux 
Comparaison des hélices 


Dans ces deux études, le moteur sera alimenté par la mise en série de trois batteries 
d'automobile de 12V, soit une tension nominale de 36V à l’entrée du contrôleur. Le choix 
de telles batteries est motivé par le fait que les batteries LiPo n'étaient pas disponibles le 
jour des essais en soufHerie. L’ensemble des essais est réalisé pour une tension aux bornes 
du contrôleur estimée à 18V (pas de mesure exacte effectuée pendant ces essais). 


Influence de la configuration de l’hélice 


Dans un premier temps, l’hélice est fixée sur le moteur de manière telle que la configu- 
ration soit tirante, la configuration usuelle d’une hélice sur un moteur (fig. 4.3, à gauche). 
Dans un second temps, l’hélice est retournée pour obtenir une configuration poussante, 
i.e. l’hélice est placée derrière le moteur (fig. 4.3, à droite). Il faudra dans ce cas inverser 
la polarité aux bornes du moteur afin d’inverser son sens de rotation. 


Le processus expérimental prévoyait de tester le moteur dans les deux configurations 
à pleine puissance. Toutefois, au vu de la quantité de courant utilisée par le moteur, 
l'élévation de température par effet Joule dans les câbles est devenue critique et certains 
contacts se sont mis à fumer. Dès lors, il a été décidé de focaliser les essais sur plusieurs 
vitesses de vent, à puissance réduite (environ 3 à 4 crans de la télécommande en deçà du 


maximum). 
. R? 
T T 
murs né 


FIGURE 4.3 — Configurations du moteur et de son hélice. 


La comparaison des courbes de poussée en fonction de la vitesse du vent pour les 
deux configurations est illustrée à la figure 4.4. Cette figure montre que la configuration 
poussante crée une force inférieure à celle produite par la configuration tirante. On peut 
expliquer cet écart par le fait que la quantité d’air aspirée dans la configuration tirante est 
légèrement supérieure à celle aspirée en configuration poussante, le moteur étant situé en 
aval de l’écoulement et non en amont comme dans le second cas. On note que l’ensemble 
des résultats suivants est obtenu avec l’hélice de 55 cm de diamètre. 
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FIGURE 4.4 - Effet de la configuration de l’hélice sur la poussée. 


Influence du diamètre de l’hélice 


Nous allons à présent étudier la différence de performance entre deux hélices iden- 
tiques mais de diamètre différent. Les deux hélices seront montées selon la configuration 
poussante, puisque c’est celle qui est privilégiée dans le design actuel du drone. Toutefois, 
la discussion précédente montre que la configuration tirante permettrait d'obtenir un gain 
de quelques newtons. 


La figure 4.5 montre la différence entre l’utilisation d’une hélice de 50 cm de diamètre et 
d’une hélice de 55 cm. Les deux hélices testées ici sont des G-Sonic du fabricant Graupner. 


Cette dernière figure montre que, comme attendu, l’hélice de moindre diamètre produit 
une poussée moindre. En effet, par conservation de la quantité de mouvement, on à que : 


= SR(Paval _ Diane) (4.1) 


où 5, est la surface de l’hélice et p sont les pressions totales en amont et en aval de l’hélice. 
Dès lors, la différence Ap étant constante, une surface d’hélice plus faible résultera dans 
une poussée moindre. 


Par vent nul, l’hélice de 55 cm permet ainsi de créer une quinzaine de newtons sup- 
plémentaires par rapport à l’hélice de 50 cm. Cependant, le gain tend à diminuer avec 
l’augmentation de la vitesse du vent. Il faut noter que nous ne disposons pas d’assez de 
mesures expérimentales pour comparer les deux hélices de manière plus détaillée. La suite 
des essais se fera naturellement avec l’hélice de 55 cm. 
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FIGURE 4.5 — Effet du diamètre de l’hélice sur la poussée. 


4.1.3 Loi de poussée 


Comme on a pu l’apercevoir au cours des paragraphes précédents, la poussée générée 
dépend de deux paramètres : la puissance fournie au moteur (action directe de l’expéri- 
mentateur sur le système) et la vitesse du vent (corollaire de la vitesse de vol désirée). 


Les essais sont cette fois réalisés avec la mise en série de deux batteries lithium- 
polymère * de 18,5V , soit une tension nominale à l'entrée du contrôleur du moteur de 
37V. Pour l’ensemble de ces tests, l’hélice et le moteur sont montés de manière à ce que 
la configuration soit poussante. 


Poussée en fonction de la puissance fournie 


Cet essai consiste à identifier la loi reliant l’énergie électrique fournie au moteur à la 
poussée délivrée par ce dernier, pour une vitesse de vent nulle. Les résultats obtenus en 
soufHerie pour deux séries de tests (avec batteries 12V et avec LiPo) sont illustrés à la 
figure 4.6. 


Cette expérience montre que l’évolution de la poussée, tant en fonction de la tension 
d’entrée que de la puissance électrique, est une courbe d'évolution quadratique. 


3. Les caractéristiques précises de ces batteries sont reprises à l’annexe D.2.2 
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FIGURE 4.6 — Poussée du moteur en fonction de la tension d’entrée ou de la puissance 
électrique. Vitesse de vent nulle. 


Poussée en fonction de la vitesse du vent 


L'intérêt de la loi de poussée en fonction de la vitesse de vent est qu’elle permet 
de déterminer les performances au décollage — nous y reviendrons à la section 4.2. En 
supposant que la poussée évolue de manière quadratique avec la vitesse, on a : 


1 = Toi + aU + œU?) (4.2) 


où 1% est la poussée statique (obtenue par vent nul), et c et « des coefficients qui sont 
déterminés par la suite. 
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FIGURE 4.7 — Poussée du moteur en fonction de la vitesse du vent, pour différentes tensions 
d'entrée. 


CHAPITRE 4 PERFORMANCES 42 


La figure 4.7 représente la poussée obtenue en soufflerie en fonction de différentes 
vitesses de vent et de plusieurs tensions d’entrée. Dès lors, partant des résultats expéri- 
mentaux de cette figure, on peut définir une loi de poussée spécifique au niveau de tension 
désiré. 


Contraintes lors des essais L’expérience montre que les batteries LiPo se déchargent 
très vite à décharge de courant maximale — c’est-à-dire si l’on impose la puissance maxi- 
male au moteur. C’est pourquoi la courbe de poussée à 24V n’est pas complète : les 
résultats obtenus en soufflerie pour des vitesses supérieures à 11 km/h étaient inaccep- 
tables et n’ont donc pas été repris. 


Calcul des coefficients c: et c2 Ces deux coefficients peuvent s’obtenir par la réso- 
lution d’un système de deux équations à deux inconnues, connaissant les coordonnées de 
deux points quelconques dans une courbe Poussée-Vitesse expérimentale. 


Intéressons-nous aux cas illustrés à la figure 4.7. Pour la courbe rouge (15 V) ainsi que 
pour la courbe verte (20 V), nous possédons assez de points afin de calculer exactement 
les deux coefficients c1 et c2. Pour le cas de la courbe bleue (24 V), nous n’avons que deux 
points, et l’un d’eux ne semble pas pertinent eu égard aux deux résultats précédents. Nous 
allons donc extrapoler les résultats, en choisissant les mêmes coefficients c, et c que ceux 
de la courbe verte, puisque cette courbe semble décroître de manière plus rapide que la 
courbe rouge (approche pessimiste). 


Les lois de poussée sont confrontées avec les essais expérimentaux dans la figure 4.8. 
On constate alors que, pour les courbes rouge et verte, la courbe de poussée obtenue 
expérimentalement et la courbe de la loi quadratique sont proches pour la plage de vitesse 
étudiée ; ceci corrobore l’utilisation d’une telle approximation quadratique. 
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FIGURE 4.8 — Loi de poussée et confrontation avec les résultats expérimentaux, pour 
différentes tensions d’entrée. 
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4.2 Performances au décollage 


L'une des phases de vol les plus critiques pour le drone est le décollage. On désire 
connaître la poussée nécessaire pour atteindre la vitesse de décollage, ainsi que l’accélé- 
ration, le temps et la longueur de piste pour y arriver. La maquette à l’échelle 1:1 du 
fuselage n’ayant pas pu être réalisée à temps, cette étude se basera sur les performances 
du moteur non installé. 


4.2.1 Poussée et accélération 


En supposant que la piste soit parfaitement horizontale et que le drone soit aligné 
parallèlement à la piste, alors la résultante des forces appliquées au drone dans l’axe de 
la propulsion, que l’on appelle T*, sera de la forme : 


T*=T — 0,59 SU°(Cn — Cr) — uW = ma (4.3) 


où Test la poussée délivrée par le moteur, a est l’accélération définie par d,U, m est la 
masse totale (constante) de l’engin et u est le coefficient de frottement inhérent au type 
de sol utilisé. 


Dans notre cas, on fera décoller le drone sur une surface dure, soit de l’asphalte ou du 
béton. La référence |[5] suggère dès lors de choisir une valeur de y égale à 0,02. 


Tous les paramètres de l’équation (4.3) dépendent de la vitesse au sol, comme par 
exemple la poussée du moteur : en effet, plus le drone accélère, plus la force de traînée 
augmente, ce qui diminue d’autant la poussée disponible. En considérant la loi de poussée 
4.2 définie précédemment, on peut calculer l’accélération instantanée à une vitesse donnée 
grâce à la relation suivante : 


mA — Jo + fiU + AU” (4.4) 


Où : 
— fo =To—-uW; 
— fi = GT; 
— fo = CTo — 0,5pS(C5 — Cr). 


L’allure de l’accélération du drone ainsi que la poussée nette (représentée pour 20 et 
24 V) sont représentées à la figure 4.9. Les paramètres utilisés pour tracer cette courbe 
sont : 

— Masse totale du drone m : 14 kg; 

— Masse volumique de l'air p : 1,225; 

— Coefficient de frottement y : 0,02; 

— Coefficient de portance Cz : 0,55 pour a — 0°; 

— Coefficient de traînée Ch : 0,07 pour a — 0°. 
Ces deux dernières valeurs sont choisies à partir des essais en soufferie reportés au chapitre 
précédent, avec une vitesse de vent de 20 m/s et un angle de calage de 2,1°. 


On peut alors observer que la traînée totale subie par le drone sur la piste est supérieure 
à la traînée en vol jusqu’à un certain point où l’inverse se produit, comme en témoigne 
une comparaison entre les figures 4.8 et 4.9(a). Ceci est dû au fait que le coefficient de 
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(a) Poussée du moteur et traînée totale. 
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(b) Accélération et force résultante selon l’axe de poussée. 


FIGURE 4.9 — Poussée, traînée et accélération du drone au décollage, en fonction de la 
vitesse d'avancement sur la piste. 


traînée utilisé dans le calcul se base sur le coefficient mesuré en soufflerie sans l'effet du 
sol. 


Par ailleurs, l’accélération sur la piste est une fonction qui dépend directement de 
l’évolution de la résultante selon l’axe de poussée : lorsque la traînée produite est égale 
à la poussée fournie par le moteur, alors l’accélération devient nulle, comme le montre la 


figure 4.9(b). 
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4.2.2 Détermination de la vitesse de décollage 


Les caractéristiques du moteur étant maintenant connues, nous pouvons déterminer 
la vitesse nécessaire au décollage. La vitesse critique au décollage est celle de décrochage ; 
ce dernier peut s’atteindre à un angle d’attaque de 6° à faible vitesse, comme l’a montré 
la figure 3.14. La vitesse de décollage VLor, définie par 1,1V,, vaut alors 13,75 m/s, et la 
vitesse de sécurité V2, définie par 1,2V,, vaut 15 m/s. 


Il est ici tentant de vérifier si la puissance délivrée par le moteur permet d’atteindre les 
vitesses proposées ci-dessus. Cependant, l’étude des performances du moteur s'effectuent 
sur le moteur seul; il s’agit donc d’une étude de la poussée non installée. Dès lors, la 
poussée mesurée expérimentalement surestime la poussée qui sera véritablement créée 
par le moteur installé sur le drone à l’échelle 1:1. Les observations suivantes sont donc 
données à titre indicatif et donnent une première idée de la poussée requise. 


En se basant sur la courbe de traînée totale au roulage tracée en noir à la figure 4.9(a), 
on remarque que la vitesse de décollage (13,75 m/s) peut être atteinte en appliquant un 
peu moins de 15 V aux bornes du moteur. 


Par ailleurs, on vérifie que la vitesse V2 peut être atteinte après le lever des roues {. 
Pour cela, on considère que le drone aura une incidence a de 5°. Dans ce cas, Cn = 0,089 
et, se basant sur la courbe de traînée totale à 5° tracée en gris à la figure 4.8, on observe 
que la vitesse V2 peut être atteinte en maintenant la tension autour de 15 V. 


Les observations précédentes montrent que tant la vitesse de décollage VLor que la 
vitesse de sécurité V2 peuvent être atteintes avec le moteur actuel sans utiliser le ré- 
gime maximum (soit plus de 24V). Cependant, ces résultats doivent être considérés avec 
précaution. Outre le fait qu'ils se basent sur des valeurs de poussée non installée, il se 
pourrait que la réserve d'énergie pour maintenir la vitesse V2 ne soit pas suffisante. Il faut 
aussi noter que la précision de la manette ne permet pas de régler finement la tension du 
moteur à distance. 


4.2.3 Longueur de piste nécessaire 


Un aspect plus pratique que l’on étudie ici est la course du drone sur la piste de 
décollage avant de décoller. On recherche donc la vitesse et la position du drone sur la 
piste en fonction du temps de roulage. Pour cela, on résout numériquement le système 


linéaire suivant : 
deyi y2 
— 4.5 
Ée 9o + g1Y2 + g2Y2 6) 


où les coefficients go, g1 et g2 correspondent respectivement aux coefficients fo, f1 et f2 
divisés par la masse totale du drone. 


Les hypothèses utilisées dans ce calcul sont : 
— le terrain est situé au niveau de la mer: 

— la vitesse du vent est nulle; 

— la température de l’air ambiant est de 15°C. 


4. S'il y en a — voir [6] 
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Les évolutions du déplacement et de la vitesse au sol du drone sont représentées à la 
figure 4.10. 
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FIGURE 4.10 — Evolutions du déplacement et de la vitesse au décollage. 


Ces graphiques montrent que, pour atteindre la vitesse VLor en appliquant 20 volts 
aux bornes du moteur, le drone devra rouler pendant moins de trois secondes, parcourant 
pendant ce temps une distance minimale de 25 mètres. Entre le début de la rotation et le 
quitter au sol du train principal, sachant que l’accélération angulaire est choisie à 11°/s 
(voir section 5.1.1) et que l’on souhaite une pente au décollage de 5 à 6°, on estime que 
la distance supplémentaire parcourue est de 10 mètres. 


En cas de problème après la vitesse de décision V1, il faut prévoir une distance suf- 
fisante ou un dispositif permettant l’arrêt du drone de manière sécurisée. La vitesse VI 
dépend de l’appareil et de son chargement. Dans le cas du drone, on supposera qu'il s’agit 
de la vitesse de décrochage V, soit 12,5 m/s; dans ce cas, le drone aura parcouru environ 
20 mètres en près de deux secondes. 


Comme le montre la figure 4.11, en l’absence de freins et en ne jouant que sur l'effet 
de la traînée, la distance parcourue jusqu’à l’arrêt complet est de plus de 180 mètres. 
Dès lors, à défaut d'utiliser une très longue piste de près de 200 mètres, il est préférable 
d'utiliser un dispositif d'arrêt, soit sur le drone lui-même (mais cette solution alourdirait 
davantage l'engin) soit dans son environnement extérieur, par exemple en installant un 
système de filet à une centaine de mètres de la ligne de départ. 


On n’oubliera pas que toutes les performances calculées ci-dessus le sont sur base 
d’une poussée non installée, et qu’elles sont donc optimistes. Il est ainsi bon de savoir 
qu’une piste de 100 mèêtres est le strict minimum au décollage, et qu’un système d’arrêt 
extérieur doit être prévu. Compte tenu de ce qui précède, on aura compris que la piste 
d'atterrissage devrait être en herbe (haute) et suffisamment longue afin de diminuer au 
maximum l'énergie cinétique avant la rencontre du filet — ou un autre dispositif à étudier. 
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FIGURE 4.11 — Distance nécessaire pour arrêter le drone après V1. 


4.3  Redéfinition du système de propulsion 


La section précédente nous a montré que la puissance actuellement disponible au dé- 
collage était adaptée, mais l’autonomie reste faible et la puissance installée reste inconnue. 
En effet, on a vu que, d’une part, la disponibilité des batteries posait rapidement pro- 
blème, sans compter les effets thermiques dus à l’ampérage important dans le circuit de 
puissance, d'autre part, la puissance réellement développée par le moteur sera moindre 
que prévue de part sa configuration poussante. 


Afin de remédier aux problèmes rencontrés ci-dessus, nous proposons deux solutions 
alternatives. La première consiste à garder la configuration monomoteur, mais en sub- 
stituant au moteur actuel un moteur plus puissant. La seconde consiste à modifier la 
configuration du système propulsif vers une configuration bimoteur. 


4.3.1 Configuration monomoteur 
Remplacement par un moteur plus puissant 


La solution la plus simple a priori semble être de remplacer le moteur actuel par un 
moteur plus puissant. Ceci permettrait de diminuer la quantité d'énergie demandée par 
le moteur, voire d'obtenir un moteur en adéquation avec le type d’hélice utilisé. Nous 
mentionnons dans la table 4.1 des moteurs brushless du catalogue Graupner susceptibles 
de convenir. Le choix de ce seul fabricant est motivé par le fait que ces moteurs sont 
reconnus dans le milieu de l’aéromodélisme comme étant les plus performants par rapport 
à leur prix. 

Cet aperçu montre que, d’une part, le choix est assez restreint — ce sont les trois seuls 
moteurs qui développent plus de 75 N nominaux dans la catégorie brushless — et d’autre 
part, le gain de puissance est tributaire d’une augmentation de poids. Le moteur 740z est 
prévu pour fonctionner avec une hélice de 55 cm et des batteries LiPo 5$, au contraire 
du moteur 6407. Il s’agit d’un choix judicieux afin d’avoir un système optimal, plutôt 
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FIGURE 4.12 - Configuration monomoteur. 


Modèle | Poussée (N) | Masse (g) | Tension nominale (V) 


6402 75 644 29,6 
740z 80 790 37 
920z 100 1708 44,4 


TABLE 4.1 — Liste de moteurs pouvant remplacer le 640z. 


qu’un moteur survitaminé avec les inconvénients qui s’ensuivent (effet Joule, autonomie 
réduite). Le moteur 9207 est quant à lui beaucoup trop lourd alors que son apport de 
puissance n’est pas nécessaire pour les vols étudiés ici. 


4.3.2 Configuration bimoteur 


Une autre solution envisageable serait de placer deux moteurs, de poussée nominale 
plus faible, sur chacun des supports sur lesquels sont fixés les poutres de l’empennage. 
Cette solution a été prévue et a été étudiée dans les travaux [2] et|6] ; elle peut donc se 
substituer à la précédente sans trop de contraintes supplémentaires. Il s’agit sans doute 
d’ailleurs de la solution la plus efficace : en effet, elle est moins tributaire de l’écoulement 
de l’air autour de la sphère, et il existe dès lors moins de risque de décollement de couche 
limite. Par ailleurs, en cas de panne moteur en vol, il existe une solution — de fortune, certes 
— de ramener le coûteux équipement au sol sans trop l’endommager. Dans la table 4.2, on 
a listé des moteurs * Graupner susceptibles de convenir à une application 2 moteurs. Les 
valeurs entre parenthèses sont les valeurs considérées en présence de 2 moteurs identiques, 
par facilité de comparaison. 


Modèle | Poussée (N) | Masse (g) | Tension nominale (V) 


555 40 (80) | 320 (640) 18,5 
555 20V | 50(100) | 320 (640) 20 
610z 48 (96) | 405 (810) 37 


TABLE 4.2 - Liste de couples de moteurs pouvant remplacer le 6407. 


Ce tableau montre immédiatement l’avantage de la nouvelle configuration proposée. 
Pour une masse quasi égale de 640 grammes, un moteur 640z développe (théoriquement) 
75 N tandis que deux moteurs 555 20V développent ensemble (théoriquement) 100 N. 


5. Le modèle 555 20V n’est plus disponible chez le fabricant en date du 1° août 2013. 
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FIGURE 4.13 - Configuration bimoteur. 


La masse des batteries n’est cependant pas prise en compte. Toutefois, si chaque moteur 
est entraîné par une seule batterie LiPo 5$S ou 6$S, de tension de sortie nominale 18,5 
ou 22,2 V, il est possible d'arriver à un bon niveau de poussée en n’augmentant pas 
substantiellement le poids du système actuel (augmentation nulle dans le cas de LiPo 5$ 
mais performances moindres, de l’ordre de 120 grammes par batterie dans le cas de LiPo 
6S mais performances maximales). 


Le choix d’une telle configuration va devoir s'accompagner d’un nouveau choix d’hélice. 
En effet, au vu de la géométrie des poutres reliant l’'empennage au fuselage (via l’aile), 
la clearance entre l’hélice et le fuselage ne permet plus d’utiliser une hélice de 55 cm 
— dans ce cas, la clearance ne serait que de 7,5 cm. De plus, une telle hélice n’est pas 
recommandée pour des moteurs de moindre puissance, puisque le couple à appliquer par le 
moteur serait trop important et diminuerait drastiquement les performances. Le fabricant 
conseille une hélice d’une quarantaine de centimètres ; une telle valeur convient également 
pour la clearance (une quinzaine de centimètres de part et d’autre du fuselage). 


La figure 4.14 montre l’allure supposée de la poussée développée par les deux moteurs 
ensemble et par un moteur seul. La loi de puissance est tracée en utilisant les mêmes 
coefficients € et c2 que ceux utilisés pour la courbe à 20 V du cas monomoteur (voir 
fig. 4.8). La courbe pour deux moteurs 555 20V, en vert, est assez proche de la courbe 
à 24 V dans le cas monomoteur. En sachant que la configuration bimoteur devrait avoir 
un meilleur rendement que la configuration monomoteur, il est supposé que la puissance 
installée sera meilleure, ce qui améliorera les performances générales. La courbe pour un 
moteur 555 20V, en bleu, montre quant à elle que la puissance devrait être suffisante pour 
maintenir le vol en palier — ce cas sera détaillé dans la section 4.4.4 sur le vol asymétrique. 
Il s’agit d’un avantage par rapport à la configuration monomoteur en cas de panne d’un 
des deux moteurs. 
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FIGURE 4.14 — Loi de poussée supposée pour 1 et 2 moteurs Graupner 555z 20V. 


4.4 Cas de vol particuliers 


Dans cette section, nous nous proposons d'étudier quatre cas de vol particuliers. Le 
premier est le cas d’une perte complète des moteurs, le second d’une perte d’un moteur 
sur deux (dans un cas bimoteur). 


4.4.1 Vol en palier 


Le vol en palier est la forme la plus simple de vol : l'appareil se déplace en ligne droite, 
à vitesse et altitude constantes. On peut alors calculer sa vitesse d’avancement optimale, 
que l’on appelle vitesse de croisière VW... Il s’agit de la vitesse à laquelle on optimise tant 
la quantité d'énergie fournie au moteur d’une part, que la traînée totale. Dans le cas 
du drone, s’agissant d’un moteur électrique, il faudra déterminer la vitesse minimale à 
laquelle le vol peut être maintenu à niveau. 


Par ailleurs, cette vitesse va dépendre de l’altitude de croisière proposée pour le vol. 
L'objectif du drone étant de prendre des prises de vue de cibles terrestres, et étant donné 
que la caméra 3D utilisée ne possède pas de zoom optique, la hauteur du vol se limitera 
à quelques dizaines de mêtres — la référence [12] conseille une hauteur de 30 mêtres. Une 
telle altitude ne va pas influencer la densité de l’air de manière apparente, donc la vitesse 
de décrochage reste 12,5 m/s et la figure 4.8 reste d'application. Cette figure montre qu’il 
est possible de voler à 20 m/s pour une tension légèrement supérieure à 15V. 


Naturellement, si l’on désire augmenter (resp. diminuer) la vitesse de vol (dans les 
limites physiques indiquées, soit V3 < U < Vxg), on devra agir sur l’élévateur afin de 
compenser le comportement cabreur (resp. piqueur) du drone. L’élévateur devra se déflé- 
chir vers le bas (resp. vers le haut) ; le pilote doit donc pousser (resp. tirer) sur le manche 
pour maintenir le vol en palier. 


Le changement des gaz devra éventuellement s'accompagner d’une correction à la 
gouverne (respectivement, vers la droite en montée et vers la gauche en plongée). Ceci 
s'applique dans le cas monomoteur, ou si les deux moteurs de la configuration bimoteur 
ne sont pas contrarotatifs. 
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4.4.2 Vol sur trajectoire circulaire 


En plus d’un maintien en palier, le drone doit pouvoir suivre une trajectoire de type 
circulaire et maintenant son altitude et sa vitesse afin de pouvoir tourner autour de la 
cible de prise de vue. On recherche donc la vitesse minimum de sustentation en virage et 
un rayon pour la trajectoire qui soit adéquat. 


Vitesse minimale en virage 
Considérons l’équilibre illustré à la figure 4.15. 
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FIGURE 4.15 - Équilibre d’un aéronef en virage. Basé sur [3]. 


En vol, la portance doit toujours équilibrer le poids de l’appareil. Lorsque l’avion fait 
un virage d’angle 6, la portance va également se déplacer d’un angle @ vers l’intérieur du 
virage. La portance de l’appareil en virage vaut L/ = L+ F,, où L est la projection de 
la portance sur l’axe vertical, tandis que F. est sa projection sur l’axe horizontal (force 
centripète). 

Pour que l’avion reste stable sur sa trajectoire, il faut donc que la force centripète soit 
compensée par une force égale et opposée (la force centrifuge F;) et que le poids soit égal 
à la composante verticale de la portance. Il suit que la portance L/ équilibre la résultante 
du poids et de la force centrifuge, appelée poids apparent : W'=W+F;. 


En conclusion, on remarque que le poids apparent est supérieur au poids propre, ce 
qui implique que la portance à créer en virage L’ doit être supérieure à la portance de vol 
en palier L. 


La portance en virage peut s’écrire en fonction de L connaissant l’angle 6 : 


L 
E = = rnb 4.6 
COS ® " 0) 


\ 5 . 1174 
où n est le facteur de charge défini par 47 — . _. 


La vitesse minimum à atteindre pour maintenir le virage s'obtient alors directement 
? 
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en se rappelant la définition du coefficient de portance et remplacant L par nW : 


2nW 
V® — _ 4.7 
: Cr. (4.7) 


La vitesse de décrochage augmente donc avec le carré du facteur de charge en virage, 
comme en témoigne la figure 4.16. Cette dernière nous apprend ainsi qu’un angle de 
virage à de plus de 50° devient dangereux, la vitesse en vol du drone étant comprise entre 
20 et 26 m/s. De plus, une telle manœuvre est difficile à maîtriser depuis le sol, puisqu'il 
faut en même temps maintenir la trajectoire, faire monter le nez de l’appareil et maintenir 
la vitesse, le tout sans afficheur d’assiette ni d’indicateur de vitesse. 


) 
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FIGURE 4.16 - Évolution de la vitesse de décrochage en fonction de l’angle du virage @. 


Rayon de la trajectoire 


On peut également déterminer le rayon de la trajectoire circulaire. On sait que l’accé- 


lération angulaire agissant sur un mobile est donnée par a = Z, soit une force centrifuge 


R ? 
Fy — Re Par ailleurs, on sait par la figure 4.15 que tan 9 = Lo ce qui donne immédia- 
tement : 


U? 
_ g tan ® 


(4.8) 


Pour une vitesse constante de 20 m/s, on peut tracer les trajectoires en fonction de 
quelques angles @, comme illustré à la figure 4.17. Cette figure nous montre que la variation 
du rayon est très rapide pour les petits angles de virage, et qu’il n’est donc pas nécessaire 
de pencher dangereusement le drone pour obtenir un rayon acceptable : un angle de 15° 
suffit pour obtenir un rayon de 150 m à 20 m/s, et ceci sans risque de s'approcher de la 
vitesse de décrochage. 
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Rayon (m) en fonction de l’angle de virage ® (°) 
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FIGURE 4.17 — Rayon de la trajectoire en fonction de quelques angles de virage @ à 20 
m/s. 


4.4.3 Pente de finesse maximale 


La méthode suivante se base sur les concepts proposés par [13]. 


Dans le cas où le système propulsif est à l’arrêt, le drone va avoir un comportement 
aérodynamique proche de celui du planeur. Il existe alors deux solutions afin de maîtriser 
sa trajectoire. La première est de rechercher des courants d’air ascendants afin de gagner 
de l’altitude au besoin, la seconde est de rechercher la pente qui permet de maximiser la 
distance franchissable avant le crash de l’appareil. C’est cette dernière solution que nous 
étudions ici. 

Soit 0 l’angle de plané. La projection du poids de l’appareil sur les axes liés à l’écou- 
lement donne : 


_— pU?SCr 
W cos 0 — (4.9) 
W sin 0 — 2 
soit, divisant les deux expressions entre elles et pour 0 petit : 
CD 
0 = — 4.10 
- (4.10) 


On peut également déterminer la vitesse du vent en fonction du poids de l’appareil à 
partir de l’expression du coefficient de portance, toujours en considérant 0 petit : 


(4.11) 
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Polaire, angle de calage (empennage) = 2,1° 
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FIGURE 4.18 — Polaire CZ vs Cn pour U — 20 m/s. 


Dans la suite, nous allons nous baser sur la polaire réalisée à une vitesse de vent de 
20 m/s et pour un angle de calage de 2,1°. Partant de l’approximation quadratique de 
cette polaire, illustrée à la figure 4.18, et connaissant l’équation de cette polaire (voir 
section 3.1.2), on peut immédiatement calculer les coefficients de traînée et de portance 
à puissance requise minimum : 


CP=:/3mX0n,. ét CAC, 


mp 


ne. ; é ; 
L’angle de plané idéal est alors donné par 0* = = ce qui permet de trouver succes- 
L 


sivement : 
— la vitesse de chute : À = —U"P6*; 
— le temps de vol restant : t,, = —£è : 
— la distance franchissable : r — — Ah. 


D 
En appliquant l’ensemble de ces relations au drone, on obtient le graphique de la dis- 
tance franchissable en fonction de la hauteur initiale, illustré à la figure 4.19. Ce graphique 
montre que la distance franchissable est bien une fonction linéaire de la hauteur de vol. 


L’angle de plané idéal du drone est aïnsi [0* — -6,6°|, ce qui donne un rapport de 
finesse f — 8,7. Une telle finesse est réaliste lorsqu'on la compare à celle d’autres avions : 

— Cessna 172 (aile haute) : f — 8,8; 

— Piper PA-28 (aile basse) : f — 10,7; 

— Jet commercial : entre 15 et 20. 
On remarque d’ailleurs dans ce dernier graphique que la distance parcourue maximale est 
près de 9 fois la hauteur de départ. 
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FIGURE 4.19 — Pente de plané maximum. 


AA.4 Vol asymétrique 


Dans le cas où le système de propulsion choisi est un système bimoteur, il importe 
d'étudier l’effet d’une panne moteur au décollage. Dans ce cas, une traînée additionnelle 
va être créée suite au frottement engendré par l’hélice immobile, ainsi qu’une plus grande 
exposition du fuselage. 


L’équation de la polaire doit être corrigée afin de prendre en compte simultanément 
l'effet de moulin à vent ainsi que la traînée supplémentaire créée par les gouvernes com- 
pensant la poussée asymétrique. 


Moulin à vent Une approximation simple de l'effet de moulin à vent peut s’écrire 
(ChSh)wm + 0,0785d3, où d2 est le diamètre d’entrée de la turbine; dans le cas d’une 
hélice, on considère que d2 = dp, le diamètre de l’hélice. 


Traînée supplémentaire Un couple de lacet va apparaître suite au déséquilibre des 
forces de part et d’autre du fuselage. Ce couple, dont le bras de levier yr vaut 50 cm, 
s'écrit : 

Ur T1 
qe Sr 
où 1; est la poussée dans le seul moteur en fonction, {F est le bras de levier de la dérive, 
q est la pression dynamique et SX est la surface totale de la dérive (dans le cas du drone, 
on aura donc deux fois la surface de la dérive). 


Nes — 


(4.12) 


La traînée totale générée par ce couple peut se diviser en deux contributions. La 
première est due à la traînée induite et est de la forme : 


MS» 
TÀF 


(CDS ind == 
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où Àr est le rapport d’allongement de la dérive. La seconde est due à la traînée de profil 
et est de la forme : 


Ajoutant cette contribution à la traînée calculée précédemment à vitesse V2, on obtient 
la polaire représentée à la figure 4.20. Cette polaire est simplement déplacée sur la droite. 


Polaire, angle de calage (empennage) = 2,1° 


0.47 —— Cas normal 
—+— Vol asymétrique 


0.05 0.1 0.15 0.2 0.25 0.3 0.35 


FIGURE 4.20 — Déplacement de la polaire vers la droite après perte moteur au décollage. 


On peut par ailleurs montrer que, malgré ce surplus de traînée et le fait qu’il ne reste 
qu’un moteur en fonction, il devrait être possible de maintenir le drone en palier, mais 
aussi d'obtenir un gradient ascendant. Ceci est représenté à la figure 4.21, établie pour 
deux angles d’attaque. 


On verra dans le chapitre suivant si la déflexion des gouvernes sur l’empennage permet 
de maintenir le nez du drone dans l’axe d’avancement souhaité. 
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en l'air 


—+— Poussée 1x 555 20V 
90F — — — Traînée totale à a = 0° 4 
— — — Traînée totale à « = 5° 


Poussée de 1 moteur et traînée totale (N) 


0 5 10 15 20 25 30 
Vitesse du vent (m/s) 


FIGURE 4.21 — Performances après lever des roues en cas de perte moteur. 


Chapitre 5 


Système de commandes 


Dans ce chapitre, nous allons nous intéresser au système de commande. Un tel sys- 
tème comprend, outre les surfaces de contrôle, l’électronique qui permet de contrôler les 
servomoteurs associés aux différentes surfaces, ainsi que le contrôleur du moteur. 


Un aspect fondamental est de choisir des surfaces de commande qui permettent de 
piloter de manière stable et précise. Par surfaces de commande ou contrôle, l’on entend 
les deux ailerons, les deux gouvernes ainsi que l’élévateur. Le drone n’est pas pourvu 
de surfaces mobiles telles que train d’atterrissage repliable, volets hypersustentateurs ou 
aérofreins. 


Un autre aspect, lié au premier, est le choix des servomoteurs à utiliser pour faire 
fonctionner les commandes. La variété de servomoteurs disponibles sur le marché est 
vaste, et un choix approprié sera proposé. 


5.1 Dimensionnement des surfaces de commande 


Les dimensions des surfaces de commande ont été étudiées dans les travaux des années 
précédentes, mais il semble opportun de réaliser une étude comparative entre les prévisions 
numériques et analytiques sur l’influence des dimensions des différentes surfaces sur le 
comportement aérodynamique du drone. 


Notons que la vitesse du drone restant relativement stable et constante au cours du 
vol, il ne semble pas opportun de le munir de dispositifs hypersustentateurs. Ceci permet 
de se passer d’une infrastructure plus complexe et d’un apport de masse non négligeable. 


5.1.1 Élévateur 


L’empennage horizontal du drone est composé d’un plan fixe et d’une commande de 
profondeur. Cette dernière, que l’on nomme élévateur, permet de contrôler le tangage de 
l’appareil, c’est-à-dire sa tendance à piquer ou à cabrer du nez. Ce type de configuration a 
été préférée à celle d’un empennage monobloc (tout l’empennage est articulé autour d’un 
axe) par sa facilité d'installation. 


Le dimensionnement de l’élévateur est tributaire du décollage. En effet, la surface doit 
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être suffisamment grande et défléchie afin d’assurer la rotation de l’appareil au sol avant 
le quitter des roues. Une comparaison avec d’autres aéronefs montrent que l’angle de 
déflexion maximum varie entre 15 et 30 degrés. 


Considérons le drone sur la piste de décollage, quelques instants avant la rotation du 
nez. L'équilibre, réalisé autour de la roue du train principal (main gear), est donné par : 


M — T(Zca + 2r) — Me + maZca + Me == Lisa (5.1) 
Où : 

— M est le moment de tangage du drone sans l’action de l’élévateur ni l’action de la 
poussée en croisière ; 

— Test la poussée nette au décollage, exprimée par l'équation (4.3), soit environ 40 
N; 

— Zca est la distance entre la piste et le centre de gravité du drone, soit 47,5 cm; 

— zr est la distance verticale entre le centre de gravité et l’axe de poussée du moteur ; 

— W est le poids du drone, soit 137 N; 

— mg est la distance selon l’axe du fuselage entre le centre de gravité du drone et son 
train principal ; 

— M, est le couple appliqué suite à la déflexion de l’élévateur ; 

— Lg est l’inertie du drone en rotation; 


— 0 est l’accélération angulaire. estimée à 11° /s. 
? 


Toutes les grandeurs précédentes sont connues, sauf le moment appliqué par l’élévateur. 
Ce dernier peut être calculé à partir de l’expression de la dérivée du coefficient de tangage 
par rapport à la déflexion de l’élévateur C»,.. Un tel coefficient peut être estimé soit de 
manière analytique, soit de manière numérique. 


Méthode analytique 


Le coefficient Cn,. peut être déterminé par un ensemble de relations analytiques ou 
empiriques issues de [14] et que l’on pourra retrouver à l’annexe C de [15]. On détermine 
par une telle méthode que pour un angle de déflexion maximal de 30°, un élévateur 
d'envergure 82 cm répond au critère de design. L’intervalle des envergures possibles est 
illustré à la figure 5.1, sur la courbe en vert. 


Méthode numérique 


Comme pour l’aérodynamique du drone, nous utilisons le code de calcul TORNADO 
pour déterminer les valeurs C,. pour plusieurs envergures de l’élévateur et deux cordes 
différentes. Les résultats sont repris à la figure 5.1. 


On remarque que, d’une part, la méthode analytique est assez fidèle au modèle numé- 
rique à corde égale. D'autre part, comme attendu, la puissance de l’élévateur s’accroît tant 
avec son envergure qu'avec sa corde. Au vu des résultats obtenus, on peut donc s’attendre 
à ce que l’élévateur remplisse convenablement sa mission en ayant les dimensions défi- 
nies actuellement, soit une envergure de 85 cm et une corde de 6 cm, pour une déflexion 
maximum de 30°. 
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—22 T T T l 

Méthode numérique, 30% corde 
—24 Méthode numérique, 40% corde 
— — — Méthode analytique, 30% corde 


_42 1 i i 


60 65 70 75 80 85 90 95 100 
Envergure de l’élévateur (cm) 


FIGURE 5.1 — Déflexion de l’élévateur en fonction des dimensions de celui-ci, vérifiant à 
la contrainte de rotation au décollage. 


5.1.2 Aiïlerons et gouvernes 


Le drone est équipé de deux dérives, chacune munie d’une gouverne, qui permettent 
de maintenir l’empennage horizontal et son élévateur en place. Les gouvernes permettent 
d'inscrire un aéronef dans son virage, ainsi que de contrer les effets de vent de travers 
pour les phases de pilotage proches du sol (décollage, atterrissage). Les gouvernes se 
défléchissent de manière symétrique. 


Par ailleurs, sur chaque aile, on trouve un aileron. Ces surfaces mobiles permettent 
de modifier la cambrure locale de l’aile et d’ainsi modifier la portance produite locale- 
ment. Les ailerons se défléchissant de manière antisymétrique, la différence de portance 
va entraîner un moment de roulis autour du centre de gravité. 


Contrairement à l’élévateur, ces deux surfaces ne peuvent se modéliser indépendam- 

ment l’une de l’autre. En effet, les effets du roulis et du lacet sont intimement liés : 

— si les gouvernes sont défléchies vers la droite, le nez de l’appareil va pointer vers la 
droite. L’extrémité de l’aile gauche va donc se déplacer plus rapidement vers la vent, 
d’où une plus grande portance, tandis que le phénomène est inversé sur l'extrémité 
de l’aile droite. Il y à alors un roulis induit, qui est favorable au virage puisqu'il 
tend à inscrire l’avion dans le virage souhaité; 

— si l’aileron de l’aile droite se défléchit vers le bas et réciproquement pour l’aile 
gauche, la quantité de portance créée sur l’aile droite est plus importante que celle 
créée sur l’aile gauche. Il y à alors apparition d’un lacet adverse, qui tend à faire 
déraper l’appareil dans la direction opposée au virage. 

Ce dernier effet peut être compensé par l'emploi d’ailerons judicieusement conçus. 
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Contraintes de design 


Les gouvernes doivent permettre un décollage ainsi qu'un atterrissage par vent de 
côté, pour un angle de lacet de 11,5° (sur proposition de [16]) ; elles permettent également 
de contrer le lacet adverse produit par les ailerons et de contrecarrer l'effet d’une perte 
moteur dans le cas d’un appareil bimoteur. La référence [17] suggère que la déflexion des 
gouvernes dans ce dernier cas doit être telle que : 


Tyr 


ME Dee, 


où yr est le bras de levier entre le moteur en fonctionnement et le centre de gravité. 


Le terme C,,, est donné par 


Ca Pie Cr VinoTe 


v 


ô 


où f/, est le coefficient volumique de la dérive, m, le rendement de la gouverne, b, l’enver- 
gure de la gouverne, b, l’envergure de la dérive et 7, l'efficacité d’angle d'attaque de la 
gouverne. Ces différentes grandeurs calculées, on obtient un angle 6, de 6,5°. 


Les surfaces mobiles doivent aussi et surtout pouvoir annuler les coefficients de moment 
de roulis et de lacet, respectivement C} et C,, soit : 
Ci = C8 + Cis, da + Cis, Or (5.2) 
Ch _ CnsP La Che, 0x + Cng, Ôr | 


Ces deux relations négligent l’effet de la traînée, de la poussée et du propwash !. Deux 
quantités parmi les précédentes sont déjà connues : les dérivées de stabilité par rapport à 
l'angle de dérapage C7, < 0 et C,, > 0. Elles sont fixées de manière à ce que la stabilité 
soit assurée. 


Résolution et discussion 


Comme pour lélévateur, les coefficients C3, Ci, Cns, et Cns, peuvent s’obtenir soit 
par des méthodes analytiques, soit numériques. Les relations analytiques sont du même 
acabit que celles utilisées pour l’élévateur. Ces coefficients peuvent aussi directement s’ob- 
tenir via un calcul statique dans TORNADO. 


Les comparaisons entre les estimations fournies par les deux méthodes pour les deux 
commandes sont illustrées à la figure 5.2 et montrent que, comme pour l’élévateur, les deux 
méthodes proposent des résultats du même ordre de grandeur. La plus grande différence se 
situe au niveau de l’annulation du lacet par l’aileron, où l’on obtient une courbe d’aspect 
parabolique jamais inférieure à 21°. 


Si l’on suit les résultats fournis par l’analyse numérique, il apparaît que les gouvernes, 
mais pas les ailerons, peuvent remplir leur mission pour un angle de déflexion d’au plus 20° 
avec les dimensions actuellement définies. Le cas des ailerons peut facilement être résolu 


1. Perturbation de la masse d’air due à l’hélice du moteur. 
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45 L 
: Lacet, 38% corde (num) 
4ot : Roulis, 42% corde (num) |] 
Lacet, 38% corde (num) 
5L a Roulis, 42% corde (num)|| 
3 ; — — — Lacet (analytique) 
: — — — Roulis (analytique) 
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Envergure d’un aileron (cm) 
(a) Déflexion d’un aïleron 
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Lacet, 30% corde (num) 
40! Roulis, 30% corde (num) |] 
dorer Lacet, 40% corde (num) 
to MEae Roulis, 40% corde (num) 
— — — Lacet (analytique) 
— — — Roulis (analytique) 


Envergure d’une gouverne (cm) 


(b) Déflexion d’une gouverne 


FIGURE 5.2 — Déflexion d’un aileron et d’une gouverne en fonction des dimensions de 
celles-ci, vérifiant le système (5.2). 


en imposant un nouvel angle limite de 25°. Cette solution est préférée à une augmentation 
de la corde de la commande, car la structure interne de l’aile doit alors être corrigée, ce qui 
n’est pas souhaitable au vu des essais déjà effectués avec fruit par [2]. Cette modification 
faite, il devient alors possible de contrôler tant le lacet que le roulis. 
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5.1.3 Dimensionnement retenu 


Pour l’élévateur, aucune modification n’est requise : l’envergure de 85 cm et la corde 
de 6 cm peuvent donc être conservées. Pour les ailerons, aucune modification n’est requise 
: l’envergure de 50 cm et la corde de 16 cm conviennent, cependant, l’amplitude du 
mouvement devrait être possible jusqu’à 25°. Pour les gouvernes, aucune modification 
n’est requise : l’envergure de 34 cm et la corde de 6 cm conviennent. 


Le dimensionnement définitif des surfaces de commande ne pourra se faire qu’à l’issue 
de tests en soufflerie sur le drone complet, afin de mesurer les dérivées des coefficients de 
moment par rapport à la déflexion des surfaces de commande. De tels essais expérimentaux 
sont encore aujourd’hui la solution utilisée par les constructeurs pour dimensionner leurs 
commandes. 


5.2 Couples appliqués aux charnières 


Les figures et équations reprises ci-après sont issues de la référence [18]. 


Nous avons étudié dans les lignes précédentes le dimensionnement des surfaces mobiles. 
Pour défléchir celles-ci, nous allons utiliser des servomoteurs de puissance et de vitesse 
appropriées. Les servomoteurs sont de petits moteurs électriques à courant continu qui 
peuvent transmettre un couple important en consommant peu d'énergie. Ils agissent di- 
rectement sur les surfaces de commande, à savoir les ailerons, les gouvernes et l’élévateur. 


Dans cette section, nous calculons les couples minimum à fournir aux charnières. Dans 
la section suivante, nous verrons comment transmettre ce couple du servomoteur à la 
commande. 


5.2.1 Calcul des dérivées des coefficients de moments aux char- 
nières 
Le couple à appliquer à une commande à s'écrit : 


U? 


où 5; est la surface de la commande et c; la corde de cette dernière. La grandeur C}, peut 
se déterminer à partir de la relation adimensionnelle : 


Ch —= ba; + b;0; (5.4) 


où à est l’angle d’attaque effectif de la commande et à l’angle de déflexion de la surface. 


Les coefficients b; et b> peuvent s'exprimer en fonction des dérivées du coefficient de 
moment par unité de longueur, comme : 


Xi 
b; . Cha — 21 x, he + AChe (5.5) 
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pour une aile sans flèche, où À; est le rapport d’allongement de la surface à et AC%,, est 
une expression définie plus loin, et 


2 
bo — Chs = Chs + ÀECha ZE, Cho + AG: (5.6) 


pour une aile sans flèche, où à, est déterminé par la figure 5.3 et AC, est une expression 
définie plus loin. 


FLAP DEFLECTION, ôg(deg) 


FIGURE 5.3 — Efficacité de portance d’une commande simple. 


Les coefficients linéiques cz, et cr, se calculent à l’aide des formules empiriques sui- 
vantes : 


Cha + th Cla g_t\lo 2)—0,5 
Cha = EE ra | Le ah tan TG a (1 Mi) (5.7) 
et 
Chs th h a g _t\|c 
… ëô € t ô ô 2\—0,5 
Chs = cn hs +2, | 1— ci (ran D +) a (1—M,) (5.8) 


Angles du bord de fuite 


Il existe plusieurs définitions de l’angle du bord de fuite. Ces définitions peuvent être 
consultées à la figure 5.4. Dans le cas qui nous occupe, on considère que la différence 
entre Dr et D est très petite, et que Pr = Pr. Dès lors, tan . vaut 0,12 pour les 
gouvernes et l’élévateur, et 0,132 pour les ailerons. 


FIGURE 5.4 — Définition des angles du bord de fuite. 
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Calcul de c;, 


Connaissant : 
— tan re ; 
— le nombre de Reynolds Re — —— : 
— le rapport cordes T ; 
— le design du nez de la commande: 
pour chaque surface portante, on peut déterminer l’ensemble des paramètres de l’équation 
(5.7) à l’aide d’un ensemble de tables et figures. Les résultats associés sont repris à la table 


9,1 


Grandeur | Aileron | Élévateur | Gouverne | Figure 
Re 6,10° J:10? 3,10° - 
Go 0,77 0,77 O,7t |C:ia) 
c® (/rad) 6,88 6,88 6,88 C.1(b) 
Cho 0,4 0,37 0,37 |C.2(a) 
ct (/rad) | -0,63 -0,53 -0,53 C.2(b) 
“No 0,8 0,8 0,8 C.3 
Ce -0,171 -0,157 -0,157 


TABLE 5.1 - Ensemble des grandeurs permettant de définir le coefficient cp... 


Calcul de b; 


Nous avons vu plus haut le terme AC... Pour une aile sans flèche, ce dernier est donné 


par la relation : 
AC 
on BK, COS Â025 


AC — C BoKa (5.9) 


où &, est la pente de la droite de portance du profil, B et le quotient sont des valeurs 
issues des graphes de la figureC.4, et K,, est une fonction de la position du contrôle sur 
le plan fixe et s'exprime par : 


LOG) RE) 
Mo — 1% 


Ke (5.10) 


où KŸ et K° sont des valeurs issues des graphes C.4, m; est le rapport de la position 
intérieure de la commande et, le même rapport pour la partie extérieure. 


Les résultats associés sont repris à la table 5.2. 


Calcul de c», 


En procédant comme précédemment, on peut déterminer l’ensemble des paramètres 
de l’équation (5.8) à l’aide d’un ensemble de tables et figures. Les résultats associés sont 
repris à la table 5.3. 
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Grandeur | Aïleron | Élévateur | Gouverne 
Ks 2.33 0,6 1 

À 7,675 5 1,85 
ACh,/X 0,005 0,008 0,018 
Ge 6,283 6,876 6,876 
B> 1 1 1 
ACha 0,073 0,033 0,123 
b; -0,063 -0,079 0,048 


TABLE 5.2 - Ensemble des grandeurs permettant de définir le coefficient b:. 


Grandeur | Aileron | Élévateur | Gouverne | Figure 
Go 0,64 0,62 0,62 C.5(a) 
Gi (/rad) Ô 4,5 4,5 C.5(b) 
Che Cha 0,76 0,78 0,78 C.6(a) 
ch, (/rad) | -0,92 -0,86 -0,86 C.6(b) 
ee a 0,7 0,7 0,7 C7 
Ce -0,46 -0,47 -0,47 


TABLE 5.3 — Ensemble des grandeurs permettant de définir le coefficient cz... 


Calcul de b 


Nous avons vu à la page 64 le terme AC. Pour une aile sans flèche, ce dernier est 
donné par la relation : 


AG 


AGE = 
Fe  BoKa COS Â025 


C5 B2K5 (5.11) 


où G,, B2 et le quotient sont des valeurs issues des graphes C.4, C.5 ou C.8, et K5 est une 
fonction de la position du contrôle sur le plan fixe et est de la même forme qu’à l’équation 


(5.10). 


Les résultats associés sont repris à la table 5.4. 


Grandeur | Aïleron | Élévateur | Gouverne 
K;5 2,33 0,6 1 
ACr,/X 0,007 0,014 0,04 
Cl, 5 46 4,5 
AChs 0,082 0,038 0,18 

b -0,4 -0,45 -0,33 


TABLE 5.4 - Ensemble des grandeurs permettant de définir le coefficient bo. 


5.2.2 Calcul de C, et H 


Comme on l’a vu à l'équation (5.4), le coefficient de moment à la charnière Cy, dépend, 
outre les coefficients b1; et b2;, des angles a; (angle d’attaque de la surface supportant la 
commande i) et 6; (angle de déflexion de la commande à). 
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L'évolution du coefficient de couple C, en fonction des angles a et à est illustrée à la 
figure 5.5, respectivement pour Ô = 0 et & = 0. Cette figure montre que l’influence de 
l’angle de déflexion de la commande a une plus grande influence sur le couple à appliquer 
que l’angle d’attaque de la surface supportant la commande. 


r r 
—— Ailerons 
—— Elévaeur 
—+- Gouvernes 


—0.05F 2 


[°2] 
{il 
Oo 


—0.1+ . 


O 
-0.15F : - 
—0.2 2 J 
nl une | 
15 —10 -5 0 5 10 15 
a (°) 
a =0 
0 T r 
—— Ailerons 
——_ Elévaeur 
—<+— Gouvernes 
—0.05 + 
—0.1+ 
T 
Cu 
-0.15F 
—0.2 
—0.25+ J 
0 5 10 15 20 25 30 


à (°) 
FIGURE 5.5 — Évolution de Cy, avec l'angle d'attaque a et l'angle de déflexion 6. 


Recherchons à présent les angles pour lesquels les contraintes sur les commandes seront 
maximales au cours du vol. 


Pour les aïlerons, l’angle a, est choisi comme étant celui donnant une portance 
minimale. Les essais réalisés en soufferie n’ont pas montré de phénomène de décrochage à 
un angle d’attaque du fuselage de -12° ; cependant, nous assumons qu'il ne saurait en être 
éloigné au vu des résultats de 2011 (-14°). L’angle d’attaque de l’aile, prenant en compte 
le calage de cette dernière, est alors choisi à -11°. La déflexion maximale de la commande 
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est, selon la mise à jour proposée à la section 5.1.3, de 25°. 


Pour l’élévateur, l’angle à. de l’empennage horizontal est, comme pour les ailerons, 
choisi pour Cr. et vaut donc, prenant en compte le calage de l’empennage (voir section 
3.4.1), -12°. La déflexion maximale de la commande est de 30°. 


Pour les gouvernes, on choisit l’angle maximal en lacet pour lequel la contrainte de 
vent latéral est vérifiée, soit 5 — 12°. La déflexion maximale des commandes est de 20°. 


Le tableau 5.5 présente les résultats obtenus pour le coefficient de moment à la char- 
nière ainsi que le couple minimal à appliquer pour permettre la déflexion des commandes 
au cours du vol. On notera que le couple est calculé à la vitesse maximale du drone, fixée 
selon la figure 4.8 à environ 26 m/s. 


Surface | Cy | H (Ncm) 
Ailerons -0,162 -86 
Élévateur | -0,22 -28 


Gouvernes | -0,105 -D,8 


TABLE 5.5 — Coefficient de moment et couple à appliquer aux charnières, selon la surface 
de contrôle. 


Ce résultat montre que la corde des commandes, plus que leur envergure, est critique 
dans le dimensionnement. Par ailleurs, si l’on ne s’autorise pas à dépasser les 20° de 
déflexion à l’aileron, alors on peut calculer que le couple maximum est de 68 Ncm, soit 
25 % du couple appliqué pour 25°. 


5.3 Transmission du couple aux commandes 


Le moment maximum aux charnières connu, il nous faut trouver un mécanisme permet- 
tant de transmettre le couple du servomoteur à la surface de commande. Les constructeurs 
proposent plusieurs solutions permettant de transmettre le couple du servo à la charnière. 
La première consiste à fixer un cylindre sur le servo qui est directement relié à la charnière, 
comme illustré à la figure 5.6(à) ; la seconde transmet le couple via un bras de levier et 
une tige rigide (poussoir) reliant le servo à une attache située sur la surface de la surface 
de commande, au niveau de la charnière. Différents types de bras sont illustrés à la figure 


5.6(b). 


5.3.1 Couples des servomoteurs 


La deuxième solution est retenue pour équiper les commandes du drone. Le couple 
calculé aux charnières sera dans ce cas différent, et dépendra du bras de levier appliqué 
tant au servomoteur, qu'entre la charnière et la surface extérieure de la commande. Soit { 
le bras de levier au niveau du servomoteur. Alors, on suppose que le couple appliqué par 
le servomoteur pour une commande donnée peut se calculer par: 


“24 
TEE 
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(a) (b) 


FIGURE 5.6 — Types de transmission pour un servomoteur. Images issues de [4] 


où L est la demi-épaisseur de la surface de commande et t’ la distance entre la surface 
et le support permettant de loger le poussoir. Il va de soi que le bras de levier £ idéal 
devrait minimiser le couple fourni par le servomoteur, tout comme le poussoir doit être de 
longueur suffisamment courte et de rigidité suffisamment élevée pour éviter un phénomène 


de flambage (voir ci-dessous). 


En utilisant un bras de levier en forme de roue et de dimensions reprises à la figure 
5.7(a), et en utilisant un support du type illustré à la figure 5.7(b) avec t’ valant environ 
1 cm, alors on peut calculer les couples appliqués par le servomoteur. Ces couples sont 
repris dans la table 5.6. 


Surface | t (mm) | { (mm) | Mervo (Ncm) 


Ailerons 19,2 10,54 -46,3 
Élévateur 7,2 10,54 7 
Gouvernes 1,2 10,54 -4,5 


TABLE 5.6 — Coefficient de moment et couple à appliquer aux servomoteurs, selon la 
surface de contrôle. 


Les servomoteurs seront choisis de manière à ce qu’ils développent le couple désiré 
pour une tension de 4,8 V, qui est la tension retenue pour l’alimentation du système de 
commandes ?. De plus, leur contrôle sera analogique. 


Après une recherche sur le site spécialisé http://wuw.servodatabase.com, on a pu 
sélectionner un ensemble de servomoteurs légers, rapides et à engrenage métallique qui 
correspondent aux couples repris ci-dessus. Cette sélection peut être consultée dans la 
table 5.7 tandis que le montage du système électronique est schématisé à la figure 5.8. 


2. Cette même tension permet d’alimenter le contrôleur du moteur. 


CHAPITRE 5. SYSTÈME DE COMMANDES 70 


FIGURE 5.7 — Bras de levier et réceptacle proposés pour connecter le servomoteur et la 
surface de contrôle au poussoir. Images issues de [1] 


Surface | Marque | Modèle | Masse (g) | Couple (Ncm) 
Aiïlerons Walkera WK-020-2 18,7 58,9 
Élévateur | XQ-Power | XQ-S1012M 12 24,5 
Gouvernes | Jamara Atom Mini 5,5 9,8 


TABLE 5.7 — Sélection de servomoteurs répondant aux critères définis. 


Aileron tribord 


Elévateur 


Gouverne tribord 


Aileron babord 


FIGURE 5.8 —- Montage du système de commandes comprenant les 5 servos et une alimen- 
tation 4,8 V. 


5.3.2 Choix des poussoirs 


On peut par ailleurs vérifier que les poussoirs, supposés idéaux, ne se déforment 


pas en flambage lorsque les forces transmises par les servomoteurs sont appliquées aux 
extrémités desdits poussoirs. Chaque surface de contrôle sera reliée à son servomoteur par 
un poussoir cylindrique plein en aluminium, de module de Young En — 70 GPa. 


La longueur des poussoirs et leur diamètre dépend de la commande. Pour les ailerons, 
on utilisera des poussoirs de 215 x @3 mm, et pour les surfaces mobiles de l’empennage, 


3. Parfaitement droits, homogènes et exempts de contrainte initiales 
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on utilisera des poussoirs de 100 x @2 mm. Ce faisant, on peut calculer la force critique 
(théorique) à partir de laquelle le phénomène de flambage se présente, et ce en utilisant 
la formule d’Euler pour une poutre rotulée à ses deux extrémités : 

T° Eaul 


Fa = — (5.12) 


A 
rd; 
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Les forces critiques calculées pour chaque commande sont reprises à la table 5.8 et 
montrent que dans tous les cas, la valeur théorique de la force menant au flambage est 
inférieure à la force appliquée par le bras du servomoteur. Il va de soi que l’utilisation de 
poussoirs non idéaux va fortement diminuer les valeurs calculées. 


où 1 — est l’inertie du poussoir et L, sa longueur. 
T 


Surface | d, (mm) | L, (mm) | Fipp (N) | Fe (N) 


Ailerons 5 215 44,1 59,5 
Élévateur 2 100 20,7 54,6 
Gouvernes 2 100 4,8 54,6 


TABLE 5.8 — Force critique de flambage des poussoirs. 


Chapitre 6 


Conclusion et Perspectives 


6.1 Conclusion 


A l'issue de ce travail, il est bon de reprendre les principaux résultats obtenus au cours 
des différentes études proposées dans les chapitres précédents. 


En ce qui concerne le fuselage, il à été expliqué l’ensemble des contraintes ayant mené 
à créer une nouvelle géométrie, et des essais en soufflerie ont été réalisés. Ces derniers ont 
démontré que les propriétés aérodynamiques étaient conservées, et même améliorées par 
rapport aux fuselages élaborés les années précédentes. On à ainsi vu que la traînée était 
moindre, et qu'il y avait une création limitée de portance. 


En ce qui concerne l’aérodynamique du drone, des comparaisons ont été effectuées entre 
différentes méthodes et des essais en soufflerie afin de capturer au mieux le comportement 
de l’appareil en vol. Suite à des problèmes techniques et par manque de temps lors des 
phases d'usinage, il n’a pas été possible de modifier le design des ailes et de l’empennage 
à tester, ni de mesurer l’ensemble des couples aérodynamiques. Cependant, on a vu que 
les forces aérodynamiques appliquées sur le modèle étaient relativement bien prédites 
par l’ensemble des méthodes testées, et que ces forces étaient du même ordre que celles 
obtenues les années précédentes. 


En ce qui concerne la propulsion et les performances du drone, nous avons montré 
que le moteur actuel était satisfaisant par rapport aux performances requises, mais qu’il 
subsistait trop d’inconnues et de défauts pour valider le concept. À nouveau, le manque 
de temps lors des phases d’usinage n’a pas permis de réaliser des essais du moteur installé. 
Des essais ont donc été réalisés sur le seul moteur afin de produire des courbes de poussée 
permettant le calcul des performances. Il à été montré que par un simple changement 
de moteur, ou par le passage à une configuration bimoteur, il était possible d'améliorer 
l'efficacité du système propulsif. L'introduction au cas bimoteur a en outre montré les 
avantages que procure une telle solution : en plus d’une puissance de propulsion plus 
élevée, il est également possible de ne voler qu'avec un moteur lors d’une panne d’un de 
ceux-ci. 


En ce qui concerne les surfaces de commande, il a été montré que les dimensions 
définies actuellement convenaient pour autant que les angles de déflexion des ailerons 
soient modifiés. Par ailleurs, les couples à appliquer aux charnières et aux servomoteurs 
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ont été calculés sur base de relations analytiques et empiriques et montrent que le couple 
le plus critique sera celui des ailerons. Enfin, le système complet de transmission du couple 
du servomoteur à la commande proposé semble répondre à nos attentes d’un point de vue 
tant mécanique que pratique. 


6.2 Perspectives 


Les points précédents montrent que ce travail n’est qu’un élément de réponse à un 
certain nombre de points encore à étudier. Parmi ceux-ci, nous proposons des pistes de 
réflexion ou des objectifs dont nos successeurs peuvent s'inspirer pour mener à bien leur 
projet. 

La première chose à faire, dès le début de l’année académique, serait de contacter M. 
Antonio Martinez afin de pouvoir construire une maquette à l’échelle 1:1 du drone. Cette 
maquette, qui peut être directement construite comme le fuselage final en fonction du 
budget, permettra de pouvoir finalement obtenir des courbes de poussée avec le moteur 
installé. Cette étude mettrait un terme à des années d’incertitude quant à la performance 
ou non du moteur actuel. 


La seconde chose, liée à la première, est de fixer définitivement le design de la confi- 
guration du (des) moteur(s). Si l’aile construite cette année est plus légère que prévu, 
il n’est pas dit qu’à la suite de la construction du fuselage, le poids soit moindre qu’es- 
compté. Il est donc intéressant, dès le début, de considérer un moteur plus puissant ou 
une configuration bimoteur. L'auteur pense que cette dernière configuration permettrait 
d'améliorer facilement l'efficacité du système propulsif tout en gagnant en redondance en 
cas de panne, le tout pour un poids égal ou légèrement supérieur au système actuel. 


En troisième étape, il serait bon d'établir les équilibres longitudinal et latéral du drone 
avec son nouveau fuselage et son empennage. Cette étude n’a pu être réalisée cette année 
suite à l’indisponibilité du capteur NANO 25. Une comparaison qualitative avec le code 
TORNADO serait également bienvenue. 


En dernier lieu, il subsiste un dernier point assez fascinant et qui nécessite des connais- 
sances approfondies en robotique ou en électronique. L'objectif initial était que le drone 
soit autonome. Il est donc souhaité qu’un système de contrôle autonome répondant à 
des critères stricts de vol et de sécurité soit développé : il faudra donc concevoir un sys- 
tème qui permet de paramétrer une altitude, une vitesse et un itinéraire. La navigation 
proposée peut se faire par programmation d’un itinéraire GPS ou par un ensemble de ba- 
lises de navigation. Le système devra aussi être pourvu d’un équipement de sécurité qui 
doit permettre de réaliser un vol exempt de tout danger pour l’environnement du drone, 
que ce soit au sol (habitations et personnes) ou en l’air (gestion du trafic aérien, TCAS, 
transpondeur, système radio, etc.) 


Annexe À 


Théorie des profils minces 


Supposons que l’aile soit infiniment longue et d'épaisseur très faible : on considère une 
plaque plane en incidence &, comme illustré à la figure A.1. 


FIGURE A.1 - Théorie des profils minces. 


On utilise une relation de transformation conforme normalisée permettant de transfor- 
mer la plaque plane en un cylindre infiniment long. Cette manipulation permet d'étudier 
l'écoulement autour du corps dans un environnement de travail connu. En effet, il est 
aisé de modéliser un cylindre par une méthode de singularité. Pour plus d'informations 
concernant une telle méthode, le lecteur pourra consulter la référence [9]. 


Une solution respectant l’imperméabilité du cylindre et la condition de Kutta! est 
d’avoir, outre un écoulement uniforme d'intensité donnée, un doublet d'intensité T7 — 
—2rU, a? ainsi qu’un tourbillon d'intensité 


F = ArU,asin(a + B) (A.1) 


où B est l’angle où se situe le bord de fuite sur le cylindre et a est le rayon du cylindre. 
Etant donné que la portance est une fonction de la circulation, on a la relation : 


E=pUT (A2) 
d’où, combinant les deux dernières équations et sachant que CC = Dr : 
ape 
8TAa . 
Cr = —sin(a + B) (A.3) 
C 


1. MA body with a sharp trailing edge which is moving through a fluid will create about itself a 
circulation of sufficient strength to hold the rear stagnation point at the trailing edge." [19] 
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Dans le cas d’une plaque mince, où € — 4a et 5 — O0, il vient immédiatement que 
Ap = °re — 27. Si cette plaque est à faible incidence, c’est-à-dire a petit, le sinus peut 


être assimilé à son angle. d’où la relation très simple pour C" : 
, L 


Cr = 27a (A.4) 


Annexe B 


Essais en soufflerie 


Cette annexe reprend un ensemble de données relatives aux essais en soufflerie. 


B.1 Calibration du mât pour la détermination des po- 
laires du nouveau fuselage 
La calibration du mât est cruciale afin de corriger la mesure et obtenir un résultat 


correct. On illustre ici la calibration utilisée lors des essais sur le drone complet à l’échelle 
1:2. 


40,00 50,00 
=|1,376x- 0,317 
35,00 - 45,00 
. y= 0,B72x+ 0,052 

40,00 

30,00 
35,00 

25,00 “ no 

20,00 25,00 

15,00 20,00 
15,00 

10,00 
10,00 

5,00 spa 

0,00 0,00 

0,00 10,00 20,00 30,00 0,00 10,00 20,00 30,00 40,00 50,00 


FIGURE B.1 — Calibration du mât dans les directions de traînée et de portance. 
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B.2 Coefficient de tangage en fonction de l’angle d’at- 
taque pour plusieurs angles de calage de l’empen- 
nage horizontal 

Lors de notre étude sur l’aérodynamique, nous nous basons sur des résultats antérieurs 
du moment de tangage, vu que le capteur NANO 25 n’était plus disponible au moment des 


tests prévus pour ce présent rapport. Les essais ont été effectués pour quatre configurations 
qui sont synthétisées dans le tableau suivant. Voir aussi [8]. 


Config. | Théorique | Mesuré 
1 -3° -1° 
2 -2° 0° 
3 d 3° 
4 3,9° 9,9" 


TABLE B.1 — Angle de calage de l’'empennage horizontal. 


C, vs «, Re -289655.1724 


0.2 
se 7 
T Config 1 
045F 2 Ne SL Config 2 
_ — Config 3 


Config 4 


K 


15 10 -5 0 5 10 
Angle d'attaque « [°] 


FIGURE B.2 - C,, en fonction de a pour plusieurs angles 44. 


B.3 Poussée pour différents régimes de vols en fonction 
de la vitesse du vent 
Résultats expérimentaux obtenus en 2011 par [15] et [20]. Ces essais ont été réalisés 


à l’aide de deux batteries LiPo 4$, et le moteur était installé sur le drone à l’échelle 1:1, 
d’où la puissance maximale moindre que celle calculée dans le présent travail. 
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70 
— Pt = 100% de Ps 
——P. = 90% deP F 
60 mot max : 
PP  =85%deP 
mot max L4 
D / 


a = 0° LA 


---D 
ma: 


x 


Poussée [N] 


V, [m/s] 


FIGURE B.3 - T'en fonction de U,, pour plusieurs cas de vols. 


Annexe C 


Calcul du couple aux charnières 


On trouvera dans les pages qui suivent les graphes utilisés pour le calcul des moments 
aux charnières du chapitre 5. Figures issues des références [14] et[18]. 


1.0 90 Y 
GRRERE nn: 
M: 9 
RER CM CY is in percent chord 


L._108 


__ SR RS 
De ER 


Note: Interpolate by plotring 
c@ 


Ce vs logjo R 
/theory N 


(tu) peory 
(per rad) 


0 04 16 20 
MAIN SUR FACETHIONESS RATIO, te 


FIGURE C.1 — Effet de l'épaisseur et de l’angle du bord de fuite sur la pente de la courbe 
de portance. 
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1.0 — 1 0—, 


(°ha) theory = 
(per rad) ° 


FIGURE C.2 — Dérivée du moment à la charnière pour une surface de commande bidimen- 


sionnelle. 
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FIGURE C.3 — Effet de la position de la charnière et de la forme du bord d’attaque sur la 
dérivée du moment. 


ERED772952D2 
HER 7474747 


[_ 
CLETITEITEITI 
ÉÉLELEDELLERE) 


FIGURE C.4 — Facteurs de correction tridimensionnels pour la dérivée du moment en 
fonction de l’angle d’attaque. 
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FIGURE C.5 — Effet de l'épaisseur et de l’angle du bord de fuite sur la pente de la courbe 
de portance. 


ANNEXE C. CALCUL DU COUPLE AUX CHARNIÈRES 83 


FIGURE C.6 — Dérivée du moment à la charnière pour une surface de commande bidimen- 
sionnelle. 
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FIGURE C.7 - Effet de la position de la charnière et de la forme du bord d’attaque sur la 
dérivée du moment. 


FIGURE C.8 — Facteurs de correction tridimensionnels pour la dérivée du moment en 
fonction de l’angle de déflexion. 


Annexe D 


Données techniques / Datasheets 


On trouvera ci-dessous l’ensemble des données techniques, datasheets et illustrations 
des appareils, composants et accessoires utilisés pour la réalisation de ce projet. 


D.1 Laboratoire des Méthodes de Fabrication 


D.1.1 Fraiseuse à commande numérique Mikron 


Ref. WF 21-D 
Course selon x 450 mm 
Course selon y 400 mm 
Course selon 2 500 mm 


Plage de vitesses | 40 — 4000 tr/min 
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D.2 Propulsion 


Toutes les données concernant la partie propulsive (moteur, contrôleur, batteries de 
puissance) sont reprises ci-dessous. 


D.2.1 Moteur brushless Graupner 6402 [1] 


Ref. 7779 
Number of poles 14 
Operating voltage range | 11,1 --- 37 V £ \ 
Nominal voltage 29,6 V 
Free shaft length 25 mm 
Case length 64 mm 
Case diameter 63 mm 
Shaft diameter 10 mm 
Weight 644 g 
No-load speed 7992 U/min 
Thrust (Propeller 50 cm) | 75 N 
Max. charging rate 55 A 
Revolutions/Volt 270 


D.2.2 Batteries LiPo RC-Plus 5$S 50C 


Ref. Python Xtreme 5000MAH 


Nombre de cellules 5 

Voltage max par cellule DU N 
Voltage min par cellule 3 V 
Courant de décharge 50C 
Courant de décharge maximum | 100C 
Charge électrique 5000 mAh 
Capacité 92,5 Wh 

| Pa PHTHOIT X TIENNE 


O2S 7.4V 570Wh "O3 11.1V 55.5Wh O4S 14.8V 740Wh 
LLLULURCPLUS.EU @ss 25.5v 25h Oess 22.2V 111.0Wh © 8S 29.SV 148.0Wh 


AIG PO E/T IL POL SIMEÈR SATTERrY 


D.2.3 Hélice Graupner G-Sonic 55 cm 


Ref. 1815.55.25 
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Diameter / Bore | 8 mm 
Diameter / Pitch | 55 x 25 cm 
Weight 200 g 


D.3 Système d’entraînement de la caméra 3D 


On trouvera ci-dessous les documents et références des moteurs, réducteurs et contrô- 
leurs suggérés pour entraîner le mécanisme de la caméra stéréoscopique. 


D.3.1 Ensemble moteur vertical : moteur 2657 W 012 CXR + 
réducteur GSR12 40:1 + contôleur IE3-512 


D 


C-Micromotors 


2 FAULHABER 
35 mNm 


Graphite Commutation 


For combination with 


Gearheads: 
26/1(S), 26A, 30/1(S), 32A 
Encoders: 
IE3-1024(L) 
Series 2657 … CXR 
2657 W 012 CXR 024 CXR 048 CXR 
1 Nominal voltage Un 12 24 48 V 
2 Terminal resistance R 0,72 2,98 12,61 Q 
3 Output power P2 max 45,3 45,7 44,1 W 
4 Efficiency, max. MN max 81 83 83 % 
5 No-load speed no 5 600 5 800 5 800 rpm 
6 No-load current (with shaft & 4 mm) lo 0,104 0,052 0,026 A 
7 Stall torque Mu 306,7 302,9 283,1 mNm 
8 Friction torque MR 2 2 2 mNm 
9 Speed constant kKn 494 247 122 rpm/V 
10 Back-EMF constant ke 2,024 4,05 8,205 mV/rpm 
11 Torque constant km 19,33 38,67 78,35 mNm/A 
12 Current constant ki 0,052 0,026 0,013 A/mNm 
13 Slope of n-M curve An/AM 18,4 19 19,6 rpm/mNm 
14 Rotor inductance L 90 365 1 500 pH 
15 Mechanical time constant Tm 3,3 3,4 3,5 ms 
16 Rotorinertia J 17 17 17 gcm? 
17 Angular acceleration © max. 180 178 172 103rad/s2 
18 Thermal resistance Rth1/Rth2 |4,4/12,6 KAW 
19 Thermal time constant Tw/Tw |28/810 s 
20 Operating temperature range: 
— motor -30 … +100 a 
- rotor, max. permissible +125 E 
21 Shaft bearings sintered bearings ball bearings, preloaded 
22 Shaft load max.: (standard) (optional version) 
— With shaft diameter 4 4 mm 
- radial at 3 000 rpm (3 mm from bearing) 10 20 N 
— axial at 3 000 rpm 2 2 N 
— axial at standstill 50 20 N 
23 Shaft play 
— radial = 0,03 0,015 mm 
— axial - 0,2 0 mm 
24 Housing material steel, zinc galvanized and passivated 
25 Weight 156 g 
26 Direction of rotation clockwise, viewed from the front face 
Recommended values - mathematically independent of each other 
27 Speed up to Ne max. 6 000 6 000 6 000 rpm 
28 Torque up to Me max. 33 34 35 mNm 
Orientation with respect to motor -0,006 0 0 -0,006 
terminals not defined 94-0010 g11 %, 05 924 $, 05 926 -0,1 910-0015 g4 -0,010 


3x 


M2 4 deep 


K2 


20,3/A 


dr 


217 


57 


| 


A © 
26,5 #7 


20,05 
0,02 


mn 7 11,6+0,5 


12403 


For notes on technical data and lifetime performance 
refer to “Technical Information”. 
Edition 2012 - 2013 


2657 W … CXR 


for Faston connector 
2,8x0,5 


© DR. FRITZ FAULHABER GMBH & CO. KG 
Specifications subject to change without notice. 


GYSIN Worm Gearbox GSR 12 


Rated torque in Nm (S1) Type Standard 
Rated torque in Nm (S5) Type Standard 
EM ERGENCY STOP torque in Nm 

Rated torque in Nm (S1) Type R Reduced 
Rated torque in Nm (S5) Type R Reduced 


EM ERGENCY STOP torque in Nm 


> 
U 
£ 
.U 
E 
ia 


” 20 2 ic 
î standard adapter plate 
Le] m) 
+ ml 
7 
. : 
Lu 
8| RIR = L = 


12 


9 19h6 


Technical and performance data: 


Gear unit weight kg 

Dimension L min. mm 

Dimension L max. mm 

Dimension | mm 

© d standard mount mm 

Lk@ max. mm 

Centre distance mm 

Average mess inertia at drive kgcm° 

No-load backlash Type Standard arcmin < 

No-load backlash Type R Reduced arcmin < 

Self-locking whenever step-down ratio > 

Max. radial load, poor fastening N 

Max. axial load N 

Torsional stiffness Nm/arcmin 

Rated speed at drive min 

Max. rated speed at drive min? 

Temperature range C -30° / +90° 
Protection class IP 44 
Service life h — 5000 


Stock programme * Standard programme + Special programme (on request) 


Encoders 


Y 


2 FAULHABER 


Magnetic Encoders 


Series IE3 — 1 024 


Lines per revolution 

Frequency range ?), up to 
Signal output, square wave 
Supply voltage 

Current consumption, typical 2) 
Output current, max. 3) 

Index Pulse width 4 

Phase shift, channel A to B 4 
Signal rise/fall time, max. (CLioap = 50 pF) 
Inertia of encoder magnet 
Operating temperature range 


D speed (rpm) = f (Hz) x60/N 
2) Uno enc = 5 V: with unloaded outputs 


Features: 

32 to 1 024 Lines per revolution 
3 Channels 

Digital output 


IE3 — 32 IE3 — 64 IE3-128 1E3-256 1E3-512 IE3 — 1024 


N 32 | 64 | 128 ES HE | 1 024 

f 15 30 60 120 240 430 kHz 
2+1 index channels 

Uoo 14,5...5,5 V DC 

ID |typ. 16, max. 23 mA 

louT |4 mA 

Po 90 + 45 EEE °e 

Co 90 + 45 90 + 75 °e 

tr/tf 10,1/0,1 us 

J 0,08 gcm? 
— 40 … + 100 *C 


3) Uno Enc = 5 V: low logic level < 0,4V, high logic level > 4,5 V: CMOS- and TTL compatible 


4) at 5 000 rpm 


Ordering information 


Encoder number lines per revolution 
of channels in combination with: 
IE3 — 32 2+1 32 DC-Micromotors 
IE3 — 64 2+1 64 2342 … CR, 2642 … CR, 2657 … CR, 
IE3 - 128 2+1 128 3242 … CR, 3257 … CR, 
IE3 - 256 2+1 256 2237 … CXR, 2642 … CXR, 2657 … CXR 
IE3 - 512 2+1 512 3272 .… CR, 3863 … CR 
IE3 — 1 024 2+1 1 024 
Brushless DC-Servomotors 
2444 … B, 3056 … B, 3564 .. B, 
4490 … B, 4490 … BS 
Features 


These incremental shaft encoders in combination with the FAULHABER 
DC-Micromotors are used for the indication and control of both shaft 
velocity and direction of rotation as well as for positioning. 


The encoder is available in a variety of different resolutions and is 
suitable for speed control and positioning applications. 
The motor and encoder are connected via separate ribbon cables. 


A permanent magnet on the shaft creates a moving magnetic field 


which is captured using a single-chip angular sensor and further 
processed. At the encoder outputs, two 90° phase-shifted rectangular 
signals are available with up to 1 024 impulses and an index impulse 


per motor revolution. 


Options 
M Connector variants (Option no.: 3807) 


AWG 28 / PVC ribbon cable (6-conductors), 
with connector PicoBlade (pitch 1,25 mm) 6 1 


M Resolutions from 1 - 127 lines per revolution 
are available by request. 


Output signals / Circuit diagram / Connector information 


Output signals 
with clockwise rotation as seen 
from the shaft end 


( P 


Amplitude 


= 


Rotation 
Admissible deviation of phase shift: 
AP= 90°- 2 x 180° £<75°  APo= |90°- 


For notes on technical data and lifetime performance 
refer to “Technical Information”. 
Edition 2012 - 2013 


Output circuit Connection Encoder 


| 


=== 1 n.c. 
2 Channell (Index) 
3 GND Enc 
4  UbbeEnc 
UDp Enc 5 ChannelB 
6 


Channel A 


GND Enc 


Caution: 

Incorrect lead connection will damage the motor 
electronics! 

When using the encoder at low temperature it 
is important to keep the cable unmoved. 


22 x 180° <75° 
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ANNEXE D. DONNÉES TECHNIQUES / DATASHEETS OL 


D.3.2 Ensemble moteur horizontal : moteur 1319 T 012SR + ré- 
ducteur 14/1 1526:1 + contrôleur IE2-400 


DC-Micromotors 


2 FAULHABER 


1,3 mNm 


Precious Metal Commutation 


Series 1319 … SR 


For combination with 
Gearheads: 

13A, 14/1, 15/5(S) 
Encoders: 

1E2-400 


024 SR 


1 Nominal voltage Un 6 12 24 
2 Terminal resistance R 8,26 34,6 119 
3 Output power P2 max 1 0,95 1,1 
4 Efficiency, max. M max 66 65 66 
5 No-load speed no 13 100 12 800 14 600 
6 No-load current (with shaft © 1,5 mm) lo 0,031 0,015 0,009 
7 Stall torque MH 2,91 2,84 2,89 
8 Friction torque MR 0,13 0,13 0,13 
9 Speed constant Kn 2 280 1110 637 
10 Back-EMF constant ke 0,438 0,897 1,57 
11 Torque constant km 4,19 8,57 15 
12 Current constant ki 0,239 0,117 0,067 
13 Slope of n-M curve An/AM 4 500 4510 5 050 
14 Rotor inductance LL 130 530 1 600 
15 Mechanical time constant in 19 19 19 
16 Rotor inertia J 0,4 0,4 0,36 
17 Angular acceleration Gta 72 71 80 
18 Thermal resistance Rth1/Rth2 |8/35 
19 Thermal time constant Tw/Tw |3,8/175 
20 Operating temperature range: 
— motor -30 … +85 (optional version -55 … +125) 
- rotor, max. permissible +125 
21 Shaft bearings sintered bearings 
22 Shaft load max.: 
— With shaft diameter 15 
- radial at 3 000 rpm (3 mm from bearing) 1,2 
— axial at 3 000 rpm 0,2 
— axial at standstill 20 
23 Shaft play 
— radial = 0,03 
— axial - 0,2 
24 Housing material steel, black coated 
25 Weight 12 
26 Direction of rotation clockwise, viewed from the front face 
Recommended values - mathematically independent of each other 
27 Speed up to Ne max. 12 000 12 000 12 000 
28 Torque up to Me max. 113 1,3 1,3 
0 0 -0,004 0 
913-0,052 96-005 g1,5-0,009 92,38 -0,04 
| A (©/20,05 |A ©) 20,07 |A 
93,5 |7|0,02 41 0,04 
a | Ï y DIN 58400 
=— m=0,2 
z—9 
x=+0,35 
1 1 
2,1 | 6+0,3 + 
lu 
19,2 |. 8,140,3 2 
1319 T … SR 1319 E .… SR 


For notes on technical data and lifetime performance 
refer to “Technical Information”. 
Edition 2012 - 2013 
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Planetary Gearheads 


2 FAULHABER 
0,3 Nm 


For combination with 
DC-Micromotors 


Series 14/1 14/1 
14/1 


Housing material 

Geartrain material 

Recommended max. input speed for: 
— continuous operation 

Backlash, at no-load 

Bearings on output shaft 

Shaft load, max.: 

- radial (6,5 mm from mounting face) 
— axial 

Shaft press fit force, max. 

Shaft play 

- radial (6,5 mm from mounting face) 
— axial 

Operating temperature range 


metal 
steel 


5 000 rpm 
<1° 
ball bearings, preloaded 


<20N 
£<5N 


Number of gear stages 1 2 3 4 5 6 
Continuous torque mNm 200 300 300 300 300 300 
Intermittent torque mNm 300 450 450 450 450 450 
Weight without motor, ca. g 17 20 24 27 30 34 
Efficiency, max. % 90 80 70 60 55 50 
Direction of rotation, drive to output = = = = = = 
Reduction ratio 3,71:1 9,7:1 43:1 94:1 415:1 2 608:1 
(rounded) 14:1 66:1 112:1 592:1 4 365:1 
134:1 989:1 5 647:1 
159:1 1 526:1 
190:1 
246:1 
L2 [mm] = length without motor 20,9 25,0 29,2 33,3 31,4 4115 
L1 [mm] = length with motor 1319T..SR 34,1 38,2 42,4 46,5 50,6 54,7 
1331T..SR 45,9 50,0 54,2 58,3 62,4 66,5 
1336U...CXR 50,9 55,0 59,2 63,3 67,4 11,5 


? The reduction ratios are rounded, the exact values are available on request or at www.faulhaber.com. 


Orientation with respect to motor 
terminals not defined 

2x 

M1,6 2,5 deep 


For notes on technical data and lifetime performance 
refer to “Technical Information”. 
Edition 2012 - 2013 


For more combinations see table. 
Example of combination with 1319...SR. 


-0,006 


0 
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TP FAULHABER 
Encoders 


Magnetic Encoders Features: 
50 to 400 Lines per revolution 
2 Channels 
Digital output 


Series IE2 — 400 


IE2 — 50 IE2 — 100 IE2 — 200 IE2 — 400 
Lines per revolution N 50 | 100 200 | 400 
Signal output, square wave 2 channels 
Supply voltage V bp 4,5... 5,5 V DC 
Current consumption, typical (Vop = 5 V DC) lon typ. 6, max. 12 mA 
Output current, max. l'our 5 mA 
Phase shift, channel À to B (o] 90 + 45 °e 
Signal rise/fall time, max. (Cioap = 50 pF) tr/tf 0,1/0,1 us 
Frequency range ?, up to f 20 | 40 | 80 | 160 kHz 
Inertia of code disc J 0,05 gcm? 
Operating temperature range — 25. + 85 ae 


9 V bp = 5 V DC: Low logic level < 0,5 V, high logic level > 4,5 V: CMOS and TTL compatible 
2 Velocity (rpm) = f (Hz) x 60/N 


Ordering information 


Encoder number lines per revolution 
of channels in combination with: 
IE2 — 50 2 50 DC-Micromotors 
IE2 — 100 2 100 1319... SR, 
IE2 — 200 2 200 1331 .… SR 
IE2 — 400 2 400 
These incremental shaft encoders in combination with the FAULHABER The supply voltage for the encoder and the DC-Micromotor as well 
DC-Micromotors are used for the indication and control of both shaft as the two channel output signals are interfaced through a ribbon 
velocity and direction of rotation as well as for positioning. cable with connector. 
The encoder is integrated in the DC-Micromotors SR-Series and Details for the DC-Micromotors and suitable reduction gearheads 
extends the overall length by only 1,7 mm! are on separate catalogue pages. 


Hybrid circuits with sensors and a low inertia magnetic disc provide 
two channels with 90° phase shift. 


Output signals / Circuit diagram / Connector information 


Output signals Output circuit Pin Function 
with clockwise rotation as seen 


from the shaft end 


[\ 
E Channel A 


1 Motor — * 
2 Motor + * 
3 GND 
p 4 Voo 
5 channel B 
6 channel A 


PVC-Ribbon cable 
6 conductors 0,09 mm? 


150 +10 


Amplitude 


œ 5/6 ———e Channel A/B 


* Note: The terminal resistance 
of all motors with precious 
metal commutation is increased 
by approx. 0.4 Q, and the max. 
allowable motor current in 
combination is 1A. Motors 
with graphite commutation 
have separate motor leads 
642 and higher motor current is 
531 allowed. 


Channel B EN 6,1 


_ 12,2 


Li 
. 
\L 


Rotation 


Admissible deviation of phase shift: 


D Connector 
A® = |90° - Pa 180° |< 45° DIN-41651 


grid 2,54 mm 


For notes on technical data and lifetime performance © DR. FRITZ FAULHABER GMBH & CO. KG 
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1319-0128SR + 14/1-1526:1 


Static - Continuous 


Series 1319-012SR + 14/1-1526:1 


1319-012SR 14/1-1526:1 


Max. voltage supply [M] 24 

Max. current supply [A] 2 

Ambient temperature [°C] 22 

Reduction [%] 0 

Motor voltage [M] 4,5 

Motor current [A] 0,051 

Motor speed [rpm] 3.052 2 
Motor torque [mNn] 0,31 260 
Motor power [W] 0,1 0,1 
Motor efficiency [%] 43 23 
Coil temperature [°C] 26,0 

Coil resistance [Q] 35,15 

Housing temperature [°C] 25,2 

Thermal resistance [K/W] 35,00 


25 
0,125 50 0,06 
_ 2,5 20 
< 0,100 E 0,05 m 
E s 2,0 d =: 
ë un Ë 0,04 3 15 a 
L d à Lu = 10 * 
S 005015 40 002 | à 
0,000 0,0 0,00 0 
0 100 200 300 400 500 600 700 800 
Torque [mNm] 
—-- Speed [rpm] - - Motor current [A] ---- Power [W] — Efficiency [%] 
Continuous limits 
3 
E 
É2 
TD 
[0] 
O 
Q 
a 1 
0 
() 50 100 150 200 250 300 
Torque [mNm] 
www.faulhaber.com 
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Annexe E 


Réglementation aérienne 


Cette annexe se propose de réunir l’ensemble des règles administratives auxquelles 
sont soumis les aéromodèles. On notera que la législation belge n’a pas, à l’heure actuelle, 
déposé de texte de loi, mais il existe une circulaire, GDF-01 4° édition de 2005, définis- 
sant un certain nombre de principes s’appliquant sur le drone que nous développons. On 
note que ladite circulaire définit également les principes légaux quant à l’utilisation et la 
propriété d’un terrain d’aéromodélisme. 


Les points issus de cette circulaire et s'appliquant à notre drone sont explicités ci-après. 


E.1 Définitions 


Un aéromodèle est un aéronef (tel que défini à l'Article 1 de la loi du 27.06.1937 
relative à la réglementation de la navigation aérienne) qui est téléguidé et utilisé à des 
fins sportives et récréatives, et qui ne transporte ni personne ni marchandise. L’aéromodèle 
est dit de catégorie 2 s’il a une masse maximale de décollage inférieure ou égale à 25 kg. 


Un terrain d’aéromodélisme est une zone bien déterminée au sol ou sur l’eau qui est 
destinée au décollage et à l’atterrissage d’aéromodèles. 


Une zone de vol est un espace aérien où l’évolution d’aéromodèles est autorisée. 


Le point de référence est le centre géométrique de la piste de décollage. Les coordon- 
nées géographiques du point de référence sont déterminées selon le Système Géodésique 
Mondial - 1984 (WGS-84 World Geodetic System, 1984) (degrés, minutes et secondes). 


E.2 Prescriptions applicables aux aéromodèles 


E.2.1 Identification 


Le drone doit être identifié par une plaque d'identification reprenant le nom et l’adresse 
complète du propriétaire. 
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E.2.2 Exigences techniques 
Niveau sonore 


Le niveau sonore émis par l’aéromodèle ne peut pas dépasser les valeurs reprises dans 
le tableau E.1. La Direction générale Transport aérien peut définir d’autres valeurs pour 
le niveau sonore compte tenu des circonstances locales existantes et compte tenu des 
législations régionales existantes. La mesure du niveau sonore est documentée en détail 
dans la circulaire. 


La distance D est la distance entre le point de référence du terrain d’aéromodélisme 
et l'habitation la plus proche. 


D (m) | Niveau sonore (dB) 


> 1200 86 
Z 900 84 
Z 600 80 


TABLE E.l -— Niveau sonore maximal autorisé. 


Mesures de sécurité 


Un aéromodéle de catégorie 2 ou supérieure doit obligatoirement être équipé d’un 
système de sécurité contre la perte de contrôle en vol suite à un brouillage électrique ou 
une panne. De plus, il doit être équipé d’un système automatique qui met les moteurs au 
ralenti ou qui coupe le moteur en cas de brouillage électrique ou de panne de l'installation 
émission et/ou réception ou en cas de brouillage des signaux radio par d’autres signaux 
radio. 


E.2.3 Procédures administratives 


Un aéromodèle de catégorie 2 ne peut être utilisé que s’il dispose d’un formulaire de 
déclaration estampillé et délivré par la DGTA. Ce formulaire de déclaration est dispo- 
nible auprès de la DGTA. L'utilisateur d’un aéromodeële de catégorie 2 doit introduire le 
formulaire de déclaration, daté et signé, en deux exemplaires à la DGTA. 


La délivrance du formulaire de déclaration estampillé est faite sur la base d’un examen 
de l’aéromodèle effectué par le responsable du terrain du club d’aéromodélisme et pour 
autant que les conditions administratives et techniques soient remplies. 


Dans le cas où une modification est apportée à l’aéromodèle, un nouveau formulaire 
de déclaration estampillé est nécessaire. 


E.3 Prescriptions applicables à la zone de vol 


La zone de vol présente une forme déterminée par les obstacles (voir circulaire pour 
plus de détail) situés aux alentours du terrain d’aéromodélisme et est limitée par un 
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cylindre dont la base est le sol, dont la hauteur maximale est de 120 m et dont le rayon 
maximal de 400 m à partir du point de référence du terrain d’aéromodélisme. Dans la 
zone de vol, les évolutions des aéromodèles doivent pouvoir être suivies continuellement à 
l’œil nu. 


E.4 Prescriptions applicables aux évolutions d’aéromo- 


dèles 


Il est interdit de faire voler des aéromodèles, entre autres: 


à partir d’un terrain pour lequel aucune autorisation d'exploitation n’a été délivrée ; 
en dehors de la zone de vol; 

si la visibilité horizontale au sol est de moins de 1000 m du point de référence du 
terrain d’aéromodélisme. La visibilité peut être confirmée par la station météo la 
plus proche ; 

à une distance de moins de 40 m du public sauf si le public se trouve derrière un filet 
de protection d’une hauteur d’au moins 2 m résistant à l’impact d’un aéromodéle ; 
au-dessus de la zone de protection du public, des constructions et des parkings ; 
au-dessus des routes secondaires quand des personnes ou des véhicules y circulent ; 
dans les nuages ; 

lorsque les évolutions des aéromodèles ne peuvent pas être suivies continuellement 
à l’œil nu ; 

à une distance de moins de 200 m d’une maison, sauf autorisation écrite et signée 
de l’occupant, et d'obstacles ; 

dans les zones de contrôle (CTR) des aéroports et dans les zones interdites (P), 
dangereuses (D) ou restrictives (R) figurant dans l’Aeronautical Information Publi- 
cation (A.[.P) sans une autorisation écrite des autorités aéronautiques qui contrôlent 
ces zones (on notera ainsi que le Sart-Tilman fait partie de la CTR de Bierset et ne 
peut donc pas être utilisé comme terrain d’aéromodélisme sans autorisation) : 

à une distance de moins de 10 km du point de référence d’un aérodrome figurant 
dans l’A.[L.P., sauf autorisation écrite de l’exploitant de l’aérodrome : 

en dehors des heures aéronautiques (entre le coucher et le lever du soleil). 
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